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ABSTRACT 

In this paper a simple analytical model for the steady-state evolution of debris disks due to collisions 
is confronted with Spitzer observations of dust around main sequence A stars. All stars are assumed 
to have planetesimal belts with a distribution of initial masses and radii. In the model disk mass 
is constant until the largest planctcsimals reach collisional equilibrium whereupon the mass falls off 
oc t~g e . Using parameters that are reasonable within the context of planet formation models and 
observations of proto-planetary disks, the detection statistics and trends seen at both 24 and 70 //m 
can be fitted well by the model. While there is no need to invoke stochastic evolution or delayed 
stirring to explain the detection statistics of dust around A stars, the model is also consistent with a 
moderate rate of stochastic events. Potentially anomalous systems are identified by their high ratio of 
observed dust luminosity to the maximum permissible in the model given their radii and ages, ///max! 
these are HD3003, HD38678, HD115892, and HD172555. It is not clear if their planetesimals have 
unusual properties (e.g., high strength or low eccentricity), or if their dust is a transient phenomenon. 
There are also well-studied examples from the literature where transient phenomena are favored (e.g., 
Vega, HD69830). However, the overall success of our model, which assumes planetesimals in all belts 
have the same strength, eccentricity and maximum size, suggests there is a large degree of uniformity 
in the outcome of planet formation. The distribution of the radii of the planetesimal belts, once 
corrected for detection bias, is found to follow N(r) cx r ~°- 8±0 - 3 i n the range 3-120 AU. Since the 
inner edge of a belt is often attributed to an unseen planet, this provides a unique constraint on 
the planetary systems of A stars. It is also shown that the effect of P-R drag on the inner edge of 
A star disks may need to be considered for those close to the Spitzer detection threshold, such as 
HD2262, HD19356, HD106591, and HD115892. Predictions are made for the upcoming SCUBA-2 
survey, including that at least 17 of the 100 A stars should be detectable above 2 mJy at 850 /xm, 
illustrating how this model can be readily applied to the interpretation of future surveys. 
Subject headings: circumstellar matter - planetary systems: formation 



1. INTRODUCTION 

The asteroid and Kuiper belts lie in regions of the so- 
lar system that are stable over long timescales to the 
gravitational perturbations of the planets (Lecar et al. 
2001). Dust produced in collisions between planetesi- 
mals in these belts inhabits a broader distribution, but 
one that is strongly influenced by the same gravitational 
perturbations (Dermott et al. 1994; Liou & Zook 1999; 
Moro-Martm & Malhotra 2003). These belts are thought 
to have started off with more than 200 times their present 

Electronic address: wyatt@ast.cam.ac.uk 

1 Michelson Science Center, California Institute of Technology, 
M/S 100-22, Pasadena, CA 91125, USA 



mass (Stern 1996; Bottke et al. 2005) and it is specu- 
lated that the belts were depleted in the same event that 
caused the Late Heavy Bombardment (LHB). The LHB 
was a period of ~ 100 Myr duration that occurred 700 
Myr after the solar system formed when the terrestrial 
planets were bombarded with a large influx of comets 
and asteroids, possibly due to restructuring of the plan- 
etary system in a dynamical instability (Gomes et al. 
2005). The evolution of these planetesimal belts since 
this event is thought to have been relatively slow, with 
mass falling off due to collisional processing, but with 
large spikes in the dust content of the inner solar system 
occurring when two large asteroids collided (Nesvorny et 
al. 2003; Farley et al. 2006). Knowledge of the plan- 
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etesimal and dust content of the solar system and how it 
evolves thus provides a rich source of information about 
both the planetesimals themselves and the planetary sys- 
tem in which they reside. 

It is important to determine the evolution of the plan- 
etesimal belts around other stars, since this may be 
equally illuminating as to the nature and evolution of 
their planetary systems. Surveys have been undertaken 
to determine this evolution, by measuring the infrared 
flux from dust produced in the collisional destruction of 
those planetesimals (Aumann et al. 1984). The observa- 
tions are used to derive basic properties of the dust belts, 
such their infrared luminosity L lT (quoted as fractional 
luminosity, / = L ir /L*), and the distance of the dust 
from the star (which is either the same as, or strongly 
dependent on, the planetesimal belt radius r) , and these 
are then compared with stellar properties, such as age 
and spectral type. Most studies have focussed on the 
evolution of dust luminosity, showing that fractional lu- 
minosity decreases with age; e.g., Spangler et al. (2001) 
found the mean dust luminosity of disks of similar age 
falls off oc i" 1 - 8 , while Greaves & Wyatt (2003) inferred 
from detection statistics that the mass falls off oc t~ - 5 , 
Najita & Williams (2005) found the mass of the detected 
disks falls off oc and the upper limit in 24 /im of the 
A star disks was found to fall off oc t~ x (Rieke et al. 
2005), while disk mass from far-IR statistics was inferred 
to be relatively age independent in the range 30-1000 
Myr (Rhee et al. 2007). These results are thus in general 
agreement with the results of theoretical models, which 
show that a planetesimal belt evolving in quasi-steady 
state would lose mass due to collisional grinding down 
giving a disk mass (and dust luminosity) that falls off 
oc r 1 (Dominik & Decin 2003; Wyatt et al. 2007). 

However, Decin et al. (2003) noted that the maximum 
luminosity of debris disks remains constant at / w 10~ 3 
even for the oldest stars, and this was explained by Do- 
minik & Decin (2003) as a consequence of delayed stir- 
ring, in which dust is not produced in the planetesimal 
belt until Pluto-sized objects form and ignite a collisional 
cascade. Since planet formation models predict that such 
massive bodies take longer to form further from the star 
(Kenyon & Bromley 2004), that interpretation makes the 
prediction that the radius of the belts should increase 
with stellar age, and more recent surveys have also con- 
sidered the evolution of radius. While in sub-mm surveys 
there is as yet no evidence of any evolution of radius with 
age (Najita & Williams 2005; Williams & Andrews 2006), 
the latest far-IR results for A stars do appear to show 
some evidence that radius is increasing with age, based 
on the fact that excess emission is detected at much later 
ages at 70 /zm than at 24 ^m (Su et al. 2006). 

It has also been suggested that much of the dust we 
see in some systems is produced episodically in collisions 
between large planetesimals, much in the same way dust 
in the zodiacal cloud increases when collisions occur in 
the asteroid belt, and that such collisions contribute to 
the large spread in the infrared luminosities seen in disks 
around A stars of similar age (Rieke et al. 2005). This 
scenario is supported by the discovery of a population of 
dust grains around Vega in the process of removal by ra- 
diation pressure which indicates that this system cannot 
have remained in steady-state for the full 350 Myr age 
of the star (Su et al. 2005). More recently it has also 



been shown that the few sun-like stars with hot dust at 
a few AU (Gaidos 1999; Beichman et al. 2005; Song et 
al. 2005; Wyatt et al. 2005) must be transient, since 
they are too bright for their age to be planetesimal belts 
that have evolved due to steady-state collisional process- 
ing; such systems may be undergoing periods analogous 
to the LHB in the solar system (Wyatt et al. 2007). 
A stochastic element to disk evolution would introduce 
additional free parameters in models and would compli- 
cate interpretation of the observed statistics if it plays a 
major role. 

It has thus not yet been possible to arrive at a con- 
sistent picture for debris disk evolution. The fact that 
surveys are undertaken for different spectral types, at dif- 
ferent wavelengths, with different sensitivities to detec- 
tions, has also added confusion to the interpretation of 
the statistics. In this paper we revisit a simple analytical 
model for the steady-state collisional evolution of plan- 
etesimal belts which was originally explored in Dominik 
& Decin (2003) and recast with slightly modified assump- 
tions in Wyatt et al. (2007). This model, including how 
it is applied in this paper to populations of debris disks, 
is given in §2. In §3 wc build a model for the population 
of disks around A stars which is constrained using the 
24 /xm statistics presented in Rieke et al. (2005). The 
properties of the model population are then compared 
in §4 with those of A star disks detected in surveys at 
70 /zm (Su et al. 2006) and the model parameters fine- 
tuned. These results indicate that, with the exception 
of a few outlier systems, the available statistics can be 
fitted by steady-state collisional evolution without hav- 
ing to invoke a major role for stochasticity or delayed 
stirring. Finally the model population is used to make 
predictions for future debris disk surveys in the sub-mm 
(§5). The conclusions are given in §6. 

2. DEBRIS DISK POPULATION MODEL 

This section describes the simple analytical model 
for the steady-state collisional evolution of planetesimal 
belts that was derived in Wyatt et al. (2007) (§2.1), and 
illustrates how that model is applied in this paper to 
populations of debris disks (§2.2). 

2.1. Steady-state collisional evolution 

The planetesimal belt is assumed to be in collisional 
equilibrium with a size distribution defined by: 

n(D) = KD 2 - 3q , (1) 

where q is 11/6 for an infinite collisional cascade 
(Dohnanyi 1969), and D is the diameter of the planetes- 
imals in km. That distribution is assumed to hold from 
the largest planetesimal in the disk, of diameter D c (in 
km), down to the size below which particles are blown 
out by radiation pressure as soon as they are created, 
Dbi (in /zm). If q is in the range 5/3 to 2 then most of 
the mass is in the largest planetesimals while the cross- 
sectional area is dominated by the smallest particles such 
that for spherical particles of density p in kg m~ 3 : 

a tot = 3.5 x 10- 17 ^(3<z - ^(W^D^f-^, (2) 

M tot = 2.5 x l0- g (3q- 5) (6 - 3q)- Vtot A,i(10 9 J D c /.Dbi) 6 i^ 

where Ctot is in AU 2 and M to t is in M ffi and the rest of 
the units are as described above and are used consistently 
throughout the paper. 



Steady-state evolution of debris disks around A stars 



3 



The planetesimal belt is assumed to be at a radius, 
r, and have a width dr, both in units of AU. Assuming 
that the grains act like black bodies, and so absorb all 
the radiation they intercept, the fractional luminosity, 
/ = Li r /L*, of the dust emission is given by 



/ = <7 to t/(47rr 2 



(4) 



In other words, in this model <7 tot , M tot , and / are all 
proportional to each other and just one is needed to de- 
fine the scaling factor K in equation (1). 

Assuming the particles act like black bodies also allows 
us to write down some other simple relations: 



D u = 0.8(L*/M*)(2700/p), 



(5) 



where L* and M* are in solar units. Also the emission 
from the disk at a wavelength A is given by: 

F„ disk = 2.35 x \Q- ll B v {\ T)a tot d~ 2 , (6) 

where d is the distance to the star in pc, giving F v in 
Jy if the Planck function B v is in Jy sr _1 . The dust 
temperature can be worked out from: 



T = 278. 3L 



0.25-0.5 



(7) 



in K. The emission from the star at the same wavelength 
is given by: 

F„ = 1.77B„(\,T+)L*T- 4 d- 2 , (8) 

where T* is in K. 

In a collisional cascade the amount of material within 
a given size range D to D + dD decreases as these plan- 
etesimals are destroyed in collisions with other members 
of the cascade, but is replaced at the same rate by frag- 
ments created from the collisional destruction of larger 
objects. The long-timescale evolution is thus determined 
by the removal of mass from the top end of the cascade. 
In this model the scaling factor K (and so the total mass 
and fractional luminosity) decreases as the number of 
planetesimals of size D c decreases. The loss rate of such 
planetesimals is determined by their collisional lifetime, 
which was derived in Wyatt et al. (2007) to be: 



3.8pr 35 (dr/r)D c \ ( {Ylq - 20) [1 + 1.25(e/J) 2 



M°-5M tot 



(18-9q)G(q,X c ) 



(9) 

in Myr, where e and I are the mean eccentricities and 
inclinations of the planetesimals' orbits, and the factor 
G(q,X c ) is given by 



G(q,X c )-. 



: [(Xt 3q - 1) + (6q - 10) (3 9 - 4)- 1 (Xf iq - 1) 

+ (3q-5)(3q-3)- 1 (Xt 3q -l)}, (10) 

and X c = D cc /D c , where D cc is the smallest planetesi- 
mal that has enough energy to catastrophically destroy 
a planetesimal of size D c . The factor X c can be worked 
out from the dispersal threshold, Q£> , defined as the spe- 
cific incident energy required to catastrophically destroy 
a particle such that (Wyatt & Dent 2002): 



X c = 1.3 x 10" 



Qt>rM; 



-I 1/3 



1.25e 2 + P 



(11) 



where Qd ^ s m J kg -1 . 

Assuming that collisions are the only process affecting 
the evolution of the disk mass M tot (f), then that mass 



(or equivalently for K, a tot, or /) can be worked out by 
solving dMtot/dt = -M to t/*c to give: 



M tot (t) = M tot (0)/[l + t/*c(0)], 



(12) 



where M to t(0) is the initial disk mass and i c (0) is the 
collisional lifetime at that initial epoch; this solution is 
valid as long as mass is the only time variable property 
of the disk. This results in a disk mass that is constant 
at Mtot(0) for t < i c (0), but that falls off oc t~ x for 
t»* c (0). 

As noted in Wyatt et al. (2007) , one property of this 
evolution is that, since the expression for t c (0) includes 
a dependence on M tot (0), the disk mass at late times is 
independent of initial disk mass. This is because more 
massive disks process their mass faster. This means that 
for any given age, t age , there is a maximum disk mass 
M max (and also infrared luminosity, /max) that can re- 
main due to collisional processing: 



M max = 



3.8 x 10" 6 pr 35 (dr/r)D 

lvl i, <-age 

(12g-20)[l + 1.25(e/J) 
(18-9q)G(q,X c ) 
10-V- 5 (dr/r) 



21-0.5 



(13) 



4ttM2-H 



ago 



2[l + 1.25(e//) 2 ] 



-0.5 



G(q,X c ) 

where again this mass is in units of M$. 



Am 
D, 



5-3? 



, (14) 



2.2. Application to populations of debris disks 

To apply this evolution to populations of debris disks, 
the following assumptions were made. All stars have a 
planetesimal belt that is undergoing a collisional cascade 
(although this is not necessarily detectable). The initial 
conditions of these belts are as described in §2.1 using 
the parameters M to t(0), r, dr, p, D Cl D Q \ (which is de- 
fined by Li, and M*) and q. The subsequent evolution 
is determined by QJ,, e and / (§2.1). To simplify the 
problem we set I — e, p = 2700 kg m -3 , dr = r/2 and 
q = 11/6 for all disks and did not consider the effects of 
changing these parameters. The parameters Qq, e and 
D c were assumed to be the same for all disks, and their 
values were constrained by a fit to the observed prop- 
erties of the known disks; the consequence of choosing 
different values is discussed. The remaining parameters, 
-^tot(O), r, L+ and M+ were chosen from distributions for 
a population of a large number of disks (10,000 in this 
case). 

For each disk the stellar spectral type was chosen at 
random from the appropriate range (i.e., B8-A9 for com- 
parison with the sample of Rieke et al. 2005), thus defin- 
ing L+, and so D a \. The distribution of initial disk 
masses, M to t(0), was based on the results of Andrews 
& Williams (2005) who did a sub-mm study of proto- 
planetary disks in the Taurus- Auriga star forming region, 
i.e., for ~ lAf© stars, from which they derived that this 
population has a log-normal distribution of dust masses 
that is centred on M mid = 3.3M e and has a la width 
of 1.14 dex. We used the same distribution but allowed 
-Mmid to be different from that found in Taurus- Auriga, 
and this was a free parameter. For the distribution of 
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disk radii, a power law distribution N(r) oc r 7 in the 
range 3-120 AU was adopted, where N(r)dr is the num- 
ber of disks with radii in the range r to r + dr. The range 
3-120 AU was chosen from the range of radii inferred for 
the sample of 46 A star disks with radius estimates (Ta- 
ble 1), and the exponent 7 was determined from a fit to 
the observations. The mass evolution of each disk was 
then completely defined using equations (9) and (12), 
and its properties at the current epoch were determined 
by assigning an age, i age , randomly from the appropri- 
ate range (i.e., 0-800 Myr for A stars). The stars were 
also assigned a random distance in the range 0-45 pc as- 
suming an isotropic distribution (i.e., N(d) oc d 2 ); given 
that there are ~ 100 A stars within 45 pc, the model is 
overpopulated by a factor of around 100 compared with 
the real population (although note that distance only be- 
comes important in this paper in §5). 

The free parameters in this population model are thus: 
M mid , 7, D c , e, and Q*. 

3. A STAR POPULATION MODEL: FIT TO 24 
STATISTICS 

Here the debris disk population model of §2 is applied 
to the statistics of the incidence of disks around A stars 
presented in Rieke et al. (2005) based on observations 
at 24 /im. There are two main aims. The first, which 
is described in the rest of this section, is to determine 
if these statistics can be reproduced with steady-state 
evolution of the disks, and if so to set constraints on the 
physical properties and starting conditions of the disks 
in that population. The second is to use the resulting 
model population to make predictions for the properties 
of A star disks found in surveys at other wavelengths and 
to determine if this simple model can explain the trends 
seen in the data. This is described in §§4 and 5, where 
the model parameters are also fine-tuned. 

3.1. Rieke et al. (2005) sample 

Rieke et al. (2005) searched a sample of 76 individ- 
ual main sequence A stars and a number of young stellar 
clusters for excess 24 /xm emission above photospheric 
levels using Spitzer and combined the data with archival 
results from IRAS for a total sample size of 266 A stars. 
They also determined the ages of all stars in a consistent 
manner thus providing a very useful sample from which 
to study the evolution of debris around A stars. Their 
main result was a plot of F 2 4 to t /-F24* versus age for all 
objects which showed a decrease in the upper envelope 
~ 150 Myr/£ age (see Fig. 1). The statistics were further 
quantified by splitting the sample into three age bins 
(< 90 Myr, 90 - 189 Myr and > 190 Myr), and classify- 
ing sources as either having small excess (i 7 24tot /Fu* < 
1.25), medium excess (1.25 < i*24tot /-F24* < 2) or large 
excess (i ? 24tot/-P , 24* > 2) (see their Fig. 3 and our middle 
Fig. 1). Data on 12 objects was also presented at 70 /im 
(see their Fig. 2 and our bottom Fig. 1 for a different 
version of this data). 

3.2. Best fit model population 

Our best fit model for the disks around a population of 
10,000 main sequence A stars is shown in Fig. 1. For each 
disk the stellar spectral type was chosen at random from 
the range B8V to A9V, with the ages of the disks chosen 
randomly from the range 0-800 Myr (that of the stars in 
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Fig. 1. — Fit of the debris disk population model to the results 
of Rieke et al. (2005) (§3). The model is comprised of 10,000 
planctcsimal belts with initial masses chosen from a log-normal 
distribution centred on lOMgj, radii in the range 3-120 AU (power 
law exponent 7 = —0.8), around stars of spectral types B8 to A9 
and ages 0-800 Myr. (Top) Total 24 fim flux divided by that of 
the stellar photosphere (cf. Fig. 1 of Rickc et al. 2005). The model 
population is shown with dots, and the disk observations of Rieke 
et al. (2005) shown with filled circles. The upper limit inferred 



from statistics in Rieke et al. (2005), p24disk/^24* 



150t a 



is shown with a solid line. (Middle) Fraction of stars with flux 
ratios in the different age bins used by Rickc et al. (2005) for 
their Fig. 3 (< 90 Myr, 90-190 Myr, > 190 Myr): F 2 4tot/-F24* = 
1-1.25 (diamond=small excess), 1.25-2 (square=medium excess), 
> 2 (triangle=large excess). Observed values are shown with y/~N 
error bars, and model values are joined with dotted, dashed and 
solid lines for increasing i<24tot /F24* ■ (Bottom) Distribution of 
radii in the model (dashed line is the whole population, dotted line 
is those predicted to be detectable at 24 and 70 fim) compared with 
the distribution inferred from observations in Rieke et al. (2005) 
with \/N error bars (asterisks). The y axis is the fraction of disks 
in each sample that fall in the different radius bins (of width 6 or 
30 AU). 
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the Rieke et al. sample) . All disks were assumed to have 
physical properties corresponding to = 300 J kg -1 , 
e = 0.05, and D c = 60 km. The distribution of their ini- 
tial masses had M m id = 10M ffi , slightly higher than that 
inferred for the disks of sun-like stars in Taurus- Auriga to 
account for a slightly higher disk mass expected around 
more massive stars (Natta et al. 2000). The consequence 
of choosing different values is discussed in §3.3. 

The distribution of radii was constrained using the 12 
disks with data at 24 and 70 /tm in Rieke et al. (2005). 
Their values of fVodisk/-F24disk were converted into dust 
temperature and then radii (both assuming black body 
emission) to derive a distribution of radii that is reason- 
ably flat, falling off only slowly between 3 and 120 AU 
(bottom Fig. 1). While this sample is small this distribu- 
tion compares well with the distribution of radii inferred 
by Najita & Williams (2005) from their sub-mm sample 
of main sequence stars, and with that inferred from the 
larger sample of 46 A star disks for which radii could be 
estimated, which is discussed in §4 (see Fig. 3). It is im- 
portant to account for the fact that this sample of stars 
with estimated radii is not representative of the whole 
population, since it only includes disks that can be de- 
tected at 24 and 70 /xm. For this reason the observed 
distribution of radii was not compared with the power 
law distribution used as input to the model, but with 
the sub-sample of disks in the model population that 
could have been detected at both 24 and 70 /im, which 
was assumed to be those with fWot /-P24* > 1-1 and 
^Votot /F70* > 1-55 (bottom Fig. 1). The exponent in the 
power law distribution was found to be 7 w —0.8 ± 0.6, 
with relatively poor constraints due to the small sample 
size. 

One immediately clear result is that the observed 
statistics can be well fitted using the model with a re- 
alistic set of planetesimal belt parameters (i.e., those 
given above) and with an initial mass distribution consis- 
tent with that inferred for proto-planetary disks. Thus 
there is no requirement for disk flux to have any stochas- 
tic component to its evolution to explain the Rieke et 
al. (2005) statistics, although this does not mean that 
stochasticity cannot play a role in determining the ob- 
servable properties of the disks. The large spread in in- 
frared excess for stars at any given age arises in the model 
from the distribution of their initial starting masses, as 
well as from that of their radii, and to a lesser extent 
from the spectral type of the parent star. Disks of high 
24 /im flux in the model are those with high initial disk 
masses and with radii that are larger than average, since 
rapid collisional processing means that close in planetes- 
imal belts tend to have low luminosities, except at the 
youngest ages. The model also predicts that disks of high 
24 /im flux should be found predominantly around higher 
luminosity stars; e.g., 40% of the stars in the model with 
^24tot /-P24* > 2 are AO or earlier compared with 25% of 
stars in the whole population. We have not been able to 
confirm this trend in the Rieke et al. (2005) data, for 
which any evidence of a change in excess ratio with spec- 
tral type is overwhelmed by the strong age dependence 
of this parameter. 

3.3. Constraints and flexibility in model parameters 

The model parameters described in §3.2 arc not 
uniquely constrained in the formulation of the evolution 



described in §2: 

(i) The upper envelope in the model population in 
top Fig. 1 is determined by / max given in equation (14). 
Thus, given the parameters and distributions that are 
fixed in the model, the upper envelope seen in top Fig. 1 
in both the model and observed populations would be 
reproduced by any combination of parameters with the 
same D^- 5 /G(q, X c ). This means that if q = 11/6 and 
X c <C 1, then the same statistics are reproduced as long 
as 

^c 1/2 Qd 5/ V 5/3 « 130 x 10 3 . (15) 

The parameters dr/r, p, e/I and q could be changed 
and accounted for by modifying the other variables in 
a similar manner. The factor derived here for A stars 
in equation (15) is four times lower than that used in 
Wyatt et al. (2007) to assess whether the hot dust of 
sun-like stars exceeds the maximum possible for steady- 
state planetesimal belts, indicating that the limit used in 
that paper is, if anything, too stringent (although it is 
possible that the properties of the planetesimal belts of 
a population of sun-like stars are different from those of 
A stars given in eq. 15). 

(ii) The mass distribution in the model determines 
the path of each disk under the upper envelope de- 
scribed above. This mass distribution thus determines 
the fraction of disks that end up in the different bins 
of i ? 24tot/-F , 24* at different times. If the relation given 
in equation (15) holds, then the same statistics as the 
model shown in Fig. 1 would also be reproduced with a 
disk mass distribution which is given by: 

M mid £ c -°- 5 « 1.3, (16) 

since, from equation (9), this would result in the same 
distribution of t c . 

In other words, if a maximum planetesimal size of 
2000 km had been assumed, the disk mass distribu- 
tion would have to be centred on M m ;d = 60M e , i.e., 
close to an order of magnitude higher than that inferred 
from the distribution of proto-planetary disks including 
non-detections (Andrews & Williams 2005), but half the 
mean mass of disks detected around Herbig Ae stars 
(Natta et al. 2000). In such an interpretation equation 
(15) indicates that e w 0.008 ^/Qp and so is satisfied by 
e = 0.1 and Q* = 140 J kg" 1 . 

Distributions with more massive disks are possible by 
assuming, e.g., that some fraction of stars end up with 
only very tenuous disks (which are not detected at 24 
/im). Lower mean disk masses are also possible to ac- 
commodate by reducing the maximum planetesimal size 
and planetesimal eccentricity; e.g., the statistics are also 
reproduced with M m ;d = IM©, D c = 0.6 km, e = 0.013, 
Qpj = 300 J kg -1 . It is not possible to differentiate 
among these possibilities at this stage. Both atypically 
high and low starting masses for debris disks have been 
suggested by other studies: e.g., a jump down in disk 
mass by a factor of 10 — 100 is inferred when extrapo- 
lating the masses of the known debris disks back in time 
and comparing with proto-planetary disk masses (Wyatt, 
Dent & Greaves 2003), and the statistics of the incidence 
of sun-like stars with planets and debris suggest that it 
is stars with disks with <~ 20./W© of solid material (i.e., 
those with masses 6 times higher than M m ;d) that end up 
with detectable debris (Greaves et al. 2006). It may also 
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be reasonable to assume that only some fraction of the 
initial disk mass ends up in the planctesimal belt, the re- 
mainder being incorporated into planets or in some other 
way removed. 

The radius distribution is also particularly poorly con- 
strained at present. On top of the large uncertainty in 
the power law exponent 7, it is also possible that the 
range of radii for the planetesimal belts is broader than 
that assumed here (3-120 AU). The model predicts that 
planetesimal belts that fall outside this range would not 
be detected in a 24 /im survey such as that of Rieke et 
al. (2005), since the rapid collisional processing of those 
< 3 AU means that these belts fall below the detection 
threshold soon after they form, while those > 120 AU 
would be too cold for their emission at 24 /im to be de- 
tectable above the photosphere. Thus we cannot rule 
out the existence of such closc-in or far-out planctesimal 
belts around up to 40% of A stars from this study. 

It is also important to point out that while the model 
presented here considered disks that all had the same 
maximum planetesimal size but with a range of starting 
masses, a model with a narrower distribution of initial 
masses and a distribution of maximum planctesimal sizes 
could also fit the observations. Likewise, the model chose 
to set planetesimal strength, eccentricity and maximum 
size to be the same for all systems. It is not a conclu- 
sion that this has to be the case, since a population for 
which these parameters are chosen from distributions can 
also fit the data. In that instance the properties of the 
planetesimal belts would also contribute to the spread in 
infrared excess seen at any given age. 

In conclusion the statistics can be readily reproduced 
with a population of disks with reasonable starting pa- 
rameters without recourse to stochastic evolution mod- 
els. The possibility remains, however, that some fraction 
of the disks in the Rickc ct al. (2005) sample are under- 
going transient events. For example, some of the disks 
may exceed the maximum luminosity permitted for disks 
at their observed radius and age. This limit is given in 
equation (14) and was used in Wyatt et al. (2007) to 
show that several of the disks found < 10 AU around 
sun- like stars must be transient. One system in the Rickc 
et al. (2005) sample that is close to this limit is £ Lep 
(HD38678), a 230 Myr A2V star which has dust with an 
infrared luminosity / = 10 x 10~ 5 (Chen & Jura 2001; 
Su et al. 2006) predicted to lie at ~ 7 AU based on its 
24 and 70 /xm fluxes, and found to have a resolved size of 
~ 3 AU at 18 /im, with emission spread across the range 
2-8 AU (Moerchen et al. 2007). For dust at 7 AU / max 
is 0.16 x 10~ 5 using the parameters inferred above for 
the rest of the A star sample (see also Table 1); i.e., this 
system is close to the limit at which we would infer it 
must be transient which was considered in Wyatt et al. 
(2007) to be when / > 100/ max . Further some fraction 
of the disks could have a two component structure, e.g., 
like the F2V star r\ Corvi which has dust at <~ 1.5 AU 
that is inferred to be transient (Wyatt et al. 2007), as 
well as at 150 AU (Wyatt et al. 2005), which is likely 
to be in steady-state. Indeed it has recently been shown 
that Vega (HD 172 167) does have an additional dust com- 
ponent at a few AU (Absil et al. 2006), as was already 
inferred for HR4796 (HD109573; Augereau et al. 1999). 
Also, there could be a stochastic element to the mid-IR 
emission which is not present in the far-IR emission, such 



as has been inferred for Vega based on the fact that the 
mid-IR emission seems to be dominated by grains under- 
going blow-out by radiation pressure (Su et al. 2005). 
Detailed study of individual objects is required to re- 
solve these uncertainties on a case-by-case basis. Since 
the threshold for detection of a stochastic event is high 
with the current observations, it is likely that the number 
of such cases is underestimated. Nonetheless, the corre- 
spondence of the time decay of the excesses with a simple 
model of steady-state evolution that has few free param- 
eters suggests that stochastic events do not dominate the 
behavior. 

3.4. Interpretation of model parameters 

Given the degeneracy in the model parameters de- 
scribed in §3.3 it is not worth dwelling on a specific in- 
terpretation of the values used in the best fit model of 
§3.2, except to say that these are reasonable parameters 
based on models of planetesimal strength (Benz & As- 
phaug 1999) and of planet formation processes (Kenyon 
& Bromley 2002). A more detailed collisional evolu- 
tion model would be required to ascertain what this is 
telling us about the physical properties of the planetesi- 
mals, for example. The important point is that the ob- 
served statistics can be explained by a model population 
in which the evolution of individual disks is very simple 
and behaves in the manner described in §2. 

4. EXPLAINING 24 AND 70 /tM A STAR RESULTS 

While there is still uncertainty in how accurately the 
simple collisional evolution model of §2.1 follows the evo- 
lution of individual disks (e.g., the validity of a single 
power law size distribution; see Wyatt et al. 2007), the 
results of §3 show that the debris disk population model 
based on that evolution can reproduce the 24 /im statis- 
tics well. In this section, the model population is further 
tested by comparing its predictions for the surveys at 
70 with the statistics observed by Su et al. (2006) 
(§4.1). The model parameters are also fine-tuned in this 
section to provide a better fit to the observed statistics 
from that paper. By plotting the model population in 
different ways we are able to provide a qualitative expla- 
nation of the trends seen in this survey, and to reproduce 
the properties of a sample of 46 A stars for which excess 
emission has been detected at both 24 and 70 /im and 
so for which estimates for the radii of their planetesi- 
mal belts can be made (§4.2). These plots provide a 
framework that can be used to interpret the statistics 
of surveys at other wavelengths, and to understand the 
implications of different detection thresholds. 

4.1. Fit to Su et al. (2006) 70 fim statistics 

The Rieke et al. (2005) 24 /im survey of A stars for 
excess emission described in §3.1 has recently been ex- 
tended to 70 /im and the results presented in Su et al. 
(2006). In that paper a sample of ~ 160 main sequence 
A-type stars (spectral types B6 to A7) were searched for 
excesses at both 24 and 70 /im using Spitzer, and ages 
were assigned to these stars using a consistent scheme. 
This sample was then supplemented with a further 19 
stars with IRAS data at 25 and 60 /im. In addition 
to reproducing the results for the evolution of 24 /im 
excesses around A stars found by Rieke et al. (2005) 
(e.g., Fig. 1), the main new result of Su et al. (2006) 
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Fig. 2. — Fit of the model population to the distributions of: (Left) 24 (im and (Right) 70 )im excesses in Su et al. (2006) (§4). The 
model parameters are unchanged to those of §3 (Fig. 1) except for a small change in planctcsimal strength. (Top) Total fluxes at 24 and 
70 /im flux divided by that of the stellar photosphere (cf. Figs. 4 and 5 of Su et al. 2006). The 10,000 disks in the model population are 
shown with dots. Stars with excess emission observed at 24 and 70 (im are shown with filled circles in the left and right plots respectively, 
while those undetected at 24 /im, or for which upper limits were derived at 70 (im, are shown with open circles. The 150ta g c une delineating 
the upper envelope of the 24 fim detections (Rickc et al. 2005) is shown with a solid line. (Bottom) Fraction of stars with flux ratios in 
the different age bins used by Su et al. (2006) for their Fig. 9 (< 30 Myr, 30-190 Myr, 190-400 Myr, > 400 Myr): F 2 Atot/F 2 A* = 1-1.25 
(diamond=small excess), 1.25-2 (square=medium excess), > 2 (trianglc=large excess); and Ffo tot / Fyo* = 1-5 (diamond=small excess), 
5-20 (square=medium excess), > 20 (triangle=large excess). Observed values are shown with y/N error bars, and model values are joined 
with dotted, dash and solid lines for increasing Ftot/F±. 



was comparable plots for the evolution of 70 /im ex- 
cesses. Plots of F to t/-P* against t ag0 for both 24 and 
70 /im (Figs. 4 and 5 of Su et al. 2006) and the frac- 
tion of stars in different age bins that have small excess 
(-p24tot/-F24* < 1-25, F 7 otot/-FVo* < 5), medium excess 
(1.25 < F 24tot /F 24 * < 2,5 < F rotot /F 70i , < 20), and 
large excess (i^tot/fk* > 2, F 70tot /F 70i , > 20) (Fig. 9 
of Su et al. 2006) are shown in Fig. 2. The observations 
show that while 70 /im excesses are also characterized by 
a t~g e fall-off, the disk emission persists much longer, and 
at much higher levels relative to the photospheric flux, 
than at 24 /im. Su et al. (2006) also identified samples 
of stars for which excess emission was detected at both 
24 and 70 fim (group I sources in Su et al.), or at 70 but 
not 24 /xm (group II sources in Su et al.), or at 24 but 
not 70 /xm (group V sources in Su ct al.). 

The model of §3 was fitted to the distribution of disk 
properties above the 24 /<m threshold, and the only way 
70 /im information was included in the model fit was 
through the 12 disks with 70 /im detections in Rieke et al. 
(2005) which set rather loose constraints on the exponent 
in the power law distribution of planetesimal belt radii 
(bottom Fig. 1). In fact it turns out that the 70 /im 
statistics of Su et al. (2006) are fitted extremely well 
using exactly the same model population. However, we 



used the new distributions of 24 and 70 /im excesses as 
a function of age shown in Fig. 2, and the distribution of 
radii for the 46 A star disks that have been detected at 24 
and 70 /im (see §4.2 and Figure 3), to set new constraints 
on the model parameters. Figures (2) and (3) show the 
resulting best fit model population, and illustrate how 
the model reproduces the longer timescale for the decay 
of 70 /im excesses seen by Su et al. (2006), and not just 
in a qualitative way: the model reproduces the observed 
statistics very well at 70 /im, just as it does at 24 /im. 
While the upper envelope in top right Fig. 2 is around 5 
times as high as that of the observed disks at any given 
age, this can be attributed to small number statistics, 
since the observed upper envelope would increase as more 
stars are included in the sample, and indeed the upper 
envelope in the model population is close to that observed 
when the number of stars included is comparable with the 
observed population (i.e., two orders of magnitude fewer 
stars in the model population than shown in Fig. 2). 

The parameters derived from this fit are almost com- 
pletely unchanged from those derived in §3, with the only 
change being in the planctcsimal strength to Qd* = 150 
J kg -1 . This means that the statistics can equally be 



8 



Wyatt et al. 




40 60 80 100 120 

Radius, AU 

Fig. 3. — Distribution of planctcsimal belt radii. That of the 
whole model population described in §4 (and shown in Fig. 2) is 
shown with a dashed line, while that of the sub-sample of the model 
population predicted to be detectable at 24 and 70 /im is shown 
with a dotted line. The observed distribution of radii for the 46 A 
star disks that have been detected at both 24 and 70 fj,m (§4.2.1) is 
shown with triangles and y/N error bars. The y axis is the fraction 
of disks in each sample that fall in the different radius bins (of 
width 6 or 30AU). 

fitted by any combination of parameters: 



£ c 1/2 Qd 5/ V 5/3 



74 x 10 3 , (17) 

which now supersedes equation (15), while the combina- 
tion of parameters required to fit the statistics that is 
given in equation (16) remains unchanged. The expo- 
nent in the radius distribution also remains unchanged 
as a result of the new fit, but this distribution is now 
better constrained to N{r) oc r~ 



-0.8±0.3 



4.2. Comparison with disks detected at 24 and 70 \xm 

In addition to the quantitative results on the incidence 
of disks of different levels of excess seen around stars of 
different ages, Su et al. (2006) identified several trends 
in their data which also require explanation in terms of 
the model. These include: (i) disks with excess at 70 /tin 
but not 24 /zm have lower fractional luminosities than 
those detected at both 24 and 70 /zm; (ii) stars that have 
excesses at 70 /im but not at 24 /zm are on average older 
than the rest of the population; 2 (iii) stars with no ap- 
parent excess have limits as low as / < 10~ 7 ; (iv) a large 
spread of two orders of magnitude of fractional luminos- 
ity at any given age; (v) older stars have lower fractional 
luminosity than younger ones; (vi) the upper envelope in 
fractional luminosity falls off oc i~ and is not constant 
as found in Decin et al. (2003); (vii) for disks detected at 
both 24 and 70 /im, the color temperature distribution 
is broader at younger ages. 

Here we explain these trends by considering the ef- 
fect of the detection threshold on the sub-samples of the 
model population that are detected at 24 or 70 //m, or 
at both wavelengths, on 3 plots of the model population. 
The model is compared with a sample of 46 A stars for 
which excess emission was detected at both 24 and 70 /im 

2 HD188228 in Su et al's group II has an age of ~ 10 Myr and 
is an exception to the rule they state that young stars with 70 (im 
excesses also have an excess at 24 fim. However, the remaining 7 
stars in their group II have ages > 325 Myr, and the mean age of 
their group II stars, ~ 490 Myr, is significantly older than the rest 
of the population. 



and so for which the radii of their disks can be estimated; 
that sample is described in §4.2.1. The three plots are 
the three fundamental disk parameters plotted against 
each other: / versus r (§4.2.2), / versus t agc (§4.2.3), 
and r versus i agc (§4.2.4). A further comparison of the 
model to the 24 and 70 /jm detected sample is given in 
§4.2.5. 

4.2.1. Sample of 46 A star disks detected at 24 and 70 

In the absence of large numbers of resolved images of 
disks from which to get a direct measure of their radii, the 
only way to estimate this important parameter is from 
interpretation of the spectrum of the excess emission, 
which must be detected at more than one wavelength to 
allow an estimate, rather than a limit, to be determined. 
Here we construct a sample of A star disks for which 
radius estimates can be made since they were all detected 
at both 24 and 70 /im. 

The sample of 46 stars is given in Table 1 and includes: 

• The 35 group I sources in Table 3 of Su et al. 
(2006) that have excesses detected at both 24 and 
70 /im (i.e., excluding HD23862 for which there is 
evidence of gas emission lines in its mid-IR spec- 
trum taken with IRS). The fractional luminosities 
of these sources were taken from Table 3 of Su et al. 
(2006), and the 24-70 /im color temperature com- 
bined with knowledge of the stellar luminosity to 
infer a planetcsimal belt radius from equation (7), 
i.e., by assuming this arises from black body emis- 
sion. Stellar ages, distances and spectral types are 
from Table 1 in Su et al. (2006), and stellar tem- 
peratures and masses were estimated from spectral 
type. 

• Eleven additional stars from the literature for 
which excess emission has been detected at both 
25 and 60 /zm by IRAS. Spitzer observations at 24 
and 70 /zm for seven of these are available either 
in the literature or in the archive, and the avail- 
able observations for these sources were analyzed 
in the same way as the Su et al. (2006) sources to 
derive stellar luminosity, dust luminosity and dust 
temperature (and so planetesimal belt radius). For 
the remaining four sources, their 25 and 60 /xm 
fluxes were taken from IRAS SCANPI 3 , and these 
fluxes were analyzed using a process analogous to 
that used in Su et al. (2006). In that process, 
spectral types were taken from SIMBAD and used 
to determine stellar temperatures and masses, and 
distances were taken from Hipparcos. Stellar lumi- 
nosities were determined by integrating under the 
spectrum of a Kurucz model atmosphere fitted to 
the K band flux from 2MASS. The stellar spec- 
trum was also used to estimate the photosphcric 
contribution to the 25 and 60 /im IRAS fluxes (hav- 
ing taken the appropriate color correction into ac- 
count) and so to determine the temperature, radius 
and fractional luminosity of the excess emission. 
Fractional luminosities and radii are consistent for 
the 3 sources that overlap with Rhee et al. (2007) 

3 http://scanpi.ipac. caltcch.edu:9000 
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TABLE 1 

Main sequence A- type stars with excess emission detected at 24 and 70 /jm (or 25 and 60 /xm) 



Star name 


Sp. Type 


T*, K M 


*, 


L* , Lq 


d, pc 


t age , Myr 


/ = L ir /L* 


r, AU 


f / /max 


///pr 


Su et al. (2006) sources 




















HD2262 


A7V 


9506 


2.9 


12.0 


23.5 


690 


0.71 x 10~ 5 


28 


1.0 


0.9 


HD 14055 


AlVnn 


9225 


2.7 


33.0 


36.1 


300 


6.7 x 10~ 5 


79 


0.65 


15 


HD19356 


B8V 


11939 


3.8 


347 


28.5 


300 


0.32 x 10 -5 


30 


1.1 


0.4 


HD27045 


A3m 


8709 


2.3 


8.0 


28.7 


193 


1.1 x 10~ 5 


22 


0.49 


1.3 


HD28355 


A7V 


7852 


1.8 


18.0 


49.2 


625 


5.6 x 10~ 5 


49 


1.8 


12 


HD30422 


A3IV 


8709 


2.3 


9.0 


57.5 


10 


4.9 x 10~ 5 


40 


0.033 


8.1 


HD31295 


A0V 


9506 


2.9 


26.0 


37.0 


10 


4.5 x 10~ 5 


78 


0.014 


9.3 


HD38056 


A0V 


9506 


2.9 


54.0 


132 


250 


4.8 x 10~ 5 


62 


0.91 


8.9 


HD38206 


A0V 


9506 


2.9 


32.0 


69.2 


9 


14 x 10~ 5 


59 


0.078 


25 


HD38678 


A2Vann 


8974 


2.5 


16.0 


21.5 


231 


9.8 x 10~ 5 


7 


100 


6.7 


HD39060 


A5V 


8203 


2.0 


13.0 


19.3 


15 


140 x 10~ 5 


24 


5.0 


200 


HD71043 


A0V 


9506 


2.9 


31.0 


73.1 


11 


4.7 x 10~ 5 


52 


0.044 


8.0 


HD71155 


A0V 


9506 


2.9 


40.0 


38.3 


169 


2.5 x 10~ 5 


44 


0.58 


3.9 


HD75416 


B8V 


11939 


3.8 


111 


96.9 


8 


4.7 x 10~ 5 


17 


0.96 


4.0 


HD79108 


A0V 


9506 


2.9 


60.0 


115 


320 


5.0 x 10~ 5 


77 


0.76 


10 


HD80950 


A0V 


9506 


2.9 


39.0 


80.8 


80 


7.5 x 10~ 5 


20 


4.7 


7.9 


HD95418 


A1V 


9225 


2.7 


68.0 


24.3 


300 


1.3 x 10~ 5 


47 


0.55 


2.1 


HD97633 


A2V 


8974 


2.5 


141 


54.5 


550 


0.71 x 10~ 5 


36 


1.4 


1.1 


HD 102647 


A3V 


8709 


2.3 


17.0 


11.1 


50 


2.0 x 10~ 5 


21 


0.40 


2.4 


HD106591 


A3V 


8709 


2.3 


26.0 


25.0 


300 


0.50 x 10~ 5 


16 


1.5 


0.5 


HD 110411 


A0V 


9506 


2.9 


23.0 


36.9 


10 


3.7 x 10~ 5 


53 


0.026 


6.3 


HD111786 


A0III 


9506 


2.9 


22.0 


60.2 


200 


3.0 x 10~ 5 


47 


0.53 


4.7 


HD 115892 


A2V 


8974 


2.5 


26.0 


18.0 


350 


1.1 x 10~ 5 


6 


32 


0.7 


HD125162 


AOp 


9506 


2.9 


16.0 


29.8 


313 


5.1 x 10~ 5 


42 


1.6 


7.7 


HD 136246 


A1V 


9225 


2.7 


31.0 


144 


15 


4.8 x 10 -5 


81 


0.021 


11 


HD 139006 


A0V 


9506 


2.9 


83.0 


22.9 


314 


1.4 x 10~ 5 


45 


0.85 


2.3 


HD158460 


A1V 


9225 


2.7 


105 


104 


260 


0.53 x 10~ 5 


68 


0.11 


1.1 


HD161868 


A0V 


9506 


2.9 


29.0 


29.1 


184 


7.5 x 10~ 5 


64 


0.72 


14 


HD165459 


A2V 


8974 


2.5 


13.0 


89.3 


5 


4.7 x 10~ 5 


31 


0.037 


6.6 


HD172167 


AOV 


9506 


2.9 


58.0 


7.8 


200 


2.3 x 10~ 5 


114 


0.093 


5.8 


HD181296 


AOVn 


9506 


2.9 


22.0 


47.7 


30 


20 x 10~ 5 


25 


2.1 


23 


HD183324 


AOV 


9506 


2.9 


22.0 


59.0 


10 


1.0 x 10~ 5 


51 


0.008 


1.8 


HD216956 


A3V 


8709 


2.3 


18.0 


7.7 


200 


6.1 x 10 -5 


67 


0.38 


13 


HD221756 


A1III 


9225 


2.7 


27.0 


71.6 


130 


1.9 X 10 


50 


0.21 


3.3 


HD225200 


AOV 


9506 


2.9 


4/.0 


129 


90 


q r\ \y in — 5 
8.0 X 10 


94 


0.20 


18 


Additional sources from literature 


















HD1438 a 


B8V 


11939 


3.8 


219 


212 


95 


ZD X 10 


114 


1.1 


/ 


HD3003 b 


B8V 


9500 


2.9 


21.0 


46.5 


50 


11 x 10~ 5 


7 


32 


6.7 


HD9672 a 


A1V 


9225 


2.7 


22.8 


61.3 


20 


72 x 10~ 5 


68 


0.54 


140 


HD21997 b 


A3IV/V 


8750 


1.9 


14.4 


73.8 


50 


47 x 10- 5 


79 


0.39 


120 


HD109573 b 


AOV 


10000 


2.9 


24.3 


67.1 


8 


330 x 10~ 5 


27 


8.7 


410 


HD110058 a 


AOV 


9506 


2.9 


10.0 


99.9 


10 


140 x 10~ 5 


18 


7.5 


140 


HD141569 b 


B9.5e 


10500 


3.1 


24.2 


99.0 


5 


460 x 10~ 5 


35 


4.4 


620 


HD153053 b 


A5IV/V 


8000 


1.8 


12.3 


50.7 


420 


7.6 x 10~ 5 


49 


1.4 


16 


HD158352 b 


A8V 


7750 


1.7 


24.0 


63.1 


600 


7.8 x 10- 5 


82 


0.85 


22 


HD172555 b 


A5IV/V 


8000 


2.0 


9.5 


29.2 


12 


51 x 10~ 5 


4 


86 


28 


HD176638 a 


B9.5V 


10000 


3.1 


47.3 


56.3 


200 


7.7 x 10~ 5 


41 


2.9 


11 



a Bascd on excess emission detected at 25 and 60 /im by IRAS 
b Based on excess emission detected at 24 and 70 fim by Spitzer 



who performed similar calculations. The ages of 
these eleven sources were taken from the literature 
(HD1438 from Wyatt et al. 2003; HD3003 from 
Song et al. 2001; HD9672, HD21997, HD110058, 
HD158352, HD172555, and HD176638 from Rhee 
et al. 2007; HD109573 from Stauffer et al. 1995; 
HD141569 from Weinberger et al. 2000; HD153053 
from Chen et al. 2006). 

For consistency we checked that the parameters given 
in Table 1 can be used to predict the observed stellar 
and excess fluxes using the equations given in this paper. 
The stellar flux should be given by equation (8), while 



the excess flux can be derived by rearranging equations 
(4), (6) and (7) to get the disk flux in Jy at a wavelength 
A to be: 



F„ disk = 2.95 x IQ- W fr 2 d- 2 B„{\,278.3L [ i- 2b r 



-0.5\ 

(18) 

Equation (18) predicts the observed disk fluxes at 24 and 
70 /zm to within 2%, which is to be expected, since these 
fluxes were used to derive the disk parameters. At 

4 The observed disk flux here is the total observed flux less the 
photospheric contribution. As such there is an uncertainty in the 
true level of disk flux due to observational uncertainties and the 
accuracy of the photospheric subtraction. 
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wavelengths other than 24 and 70 fxm, there may be in- 
accuracies due to the assumption that the emission spec- 
trum is described by that of a black body at a single 
temperature (e.g., see §5). The resulting stellar fluxes 
are also in general agreement, although they are for some 
stars in error by a factor of up to 2. We attribute this 
to the shape of the stellar spectrum, which means that a 
bolometric luminosity cannot be rigidly applied to derive 
stellar fluxes to a greater accuracy, but consider that this 
accuracy is sufficient to perform a simple analysis of this 
sample. 

Since some of these disks have been resolved in ground- 
and space-based imaging, it is worth comparing the re- 
solved sizes with those predicted by this crude analysis 
of the 24 and 70 /im fluxes. For the disks seen to be con- 
fined to rings it is clear that the prediction given in Table 
1 underestimates the true size by a factor of 2 — 3: Foma- 
lhaut is predicted to have dust at 67 AU whereas 133 AU 
is its observed radius (Kalas et al. 2005), and HR4796 
is predicted to have dust at 27 AU whereas 70 AU is 
observed (Telesco et al. 2000). Likewise HD181327 is re- 
ported to have an observed size that is 3 times larger than 
that predicted using black body grains (Schneider et al. 
2006). This discrepancy can be readily explained by the 
fact that the smallest dust grains in their distributions, 
which dominate the cross-sectional area, emit at a higher 
temperature than black body for a given distance from 
the star (Wyatt et al. 1999; Wyatt & Dent 2002). For the 
two young systems HD39060 ((3 Pictoris) and HD141569, 
the dust is observed to span a large range of radii: from 
tens of AU out to 1000s of AU for (3 Pictoris with a den- 
sity distribution that peaks at 75 AU (Kalas & Jewitt 
1995; Telesco et al. 2005), and from - 10 to 1200 AU 
for HD141569 (Fisher et al. 2000; Clampin et al. 2003). 
These extended distributions certainly complicate the in- 
terpretation of these systems in terms of rings, although 
it is believed that radiation forces cause the distributions 
of dust in these systems to extend far beyond those of 
the planetesimal belts that created them (Augereau et 
al. 2001), so that the predicted planetesimal belt radii 
of 24 AU for /? Pictoris and 35 AU for HD141569 arc not 
as inaccurate as they first appear. For Vega, the pre- 
dicted size of 114 AU is close to the size observed in the 
sub- mm (Holland et al. 1998), even if the far-infrared 
emission on which the size estimate is based extends out 
to much larger distances (Su et al. 2005). We have opted 
to retain the radius derived from the SED modeling for 
plotting purposes, since these are then derived in a more 
consistent manner across the dataset. 

All but four of the disks in this sample were detected by 
Spitzer at both 24 and 70 /im. The remaining four disks 
have yet to be observed by Spitzer, but were detected by 
IRAS at both 25 and 60 fim. Thus we consider that the 
sample derived here is representative of those disks which 
it is possible to detect at 24 and 70 ^m, and in the follow- 
ing sections we compare this sample with the sub-sample 
of the model population that would have been detected in 
the study of Su et al. (2006), assuming the detection lim- 
its quoted by those authors. The limiting factor for these 
observations was not instrumental sensitivity, rather it 
was the ability to distinguish between photospheric and 
excess emission. It was thus the accuracy of the extrap- 
olation of the photospheric flux from near-IR to longer 
wavelengths that determined the detection threshold. A 



detection required a disk flux relative to the stellar flux 
at 24 and 70 ^m of ^disk/i^* > i?24dot = 0.1 (for a 
5cr detection) and FVodisk/^Vo* > -fiVodet = 0.55 (for a 3cr 
detection), respectively. Rearranging equations (8) and 
(18) shows that such detection limits correspond to disks 
with a fractional luminosity of: 

/dot = 6.0xl0 9 i? dct r- 2 ^T- 4 B l ,(A,r,)/B 1/ (A,278.3^- 25 r- - 5 ). 

(19) 

4.2.2. Fractional luminosity vs radius 

The most valuable plot for interpreting the statistics 
of debris disk surveys that have detection thresholds 
given by equation (19) is the plot of fractional luminos- 
ity against radius, Fig. 4. The 10,000 disks in the model 
population are shown in left Fig. 4. There is a general 
trend evident in the upper envelope of the model popu- 
lation following / oc r 7 / 3 . This trend can be understood 
from equation (14), and considering that when X c <C 1 
then G(ll/6, X c ) oc r~ 5 / 6 , and so / max oc r 7 / 3 (Wyatt et 
al. 2007). Lines showing the maximum possible luminos- 
ity for planetesimal belts in the model that have ages of 
8, 80 and 800 Myr, assuming these are around A0V stars 
with the parameters for QJ,, e and D c described in §3, are 
also shown on this figure. These lines are approximately: 

/ max « 1.2 x 10-V/ 3 ^. (20) 

In other words, the maximum possible luminosity for 
planetesimal belts of a given radius depends only on age, 
and the upper envelope in the model population is thus 
determined by the youngest age of its disks. 

The lines of maximum luminosity are extremely useful, 
since the evolution on Fig. 4 of any one planetesimal belt 
in the model, with initial parameters of radius r\ and 
fractional luminosity /1, is for it to remain stationary at 
(n, /1) while its age is i ago < 1.2 x 10" 6 r 7 / 3 /f 1 Myr, 
and then for its fractional luminosity to decrease at later 
times such that it is always found at (n, /max) at an 
age i ago > 1.2 x lO-V/Vf 1 Myr. This means that a 
disk of a given age either lies on the appropriate line 
of maximum luminosity, or below it, so that, e.g., the 
sub-sample of model disks that are older than 80 Myr 
are (some fraction of) the disks that lie below the line 
of maximum luminosity for 80 Myr. Note that while 
this is useful for illustrating the trends in the data it 
should not be applied too rigidly, since these lines are 
also a function of spectral type, and their exact location 
depends on parameters like QJ,, e and D c , which may 
vary among disks. 

Also shown on left Fig. 4 is the detection limit at 24 
and 70 fim appropriate for planetesimal belts around 
A0V stars assuming the detection thresholds given in 
§4.2.1 (equation 19). The detection threshold at 24 /im 
is roughly flat at /24dot = 0.5 — 2 x 10~ 5 across all radii, 
while that at 70 /im falls more than an order of magni- 
tude across the range 3-120 AU. Disks in the model pop- 
ulation that we predict could be detected at 24 and 70 
/im would lie above both of these lines (noting, however, 
that the location of the line is also spectral type depen- 
dent). Such a sub-sample is shown in right Fig. 4 where 
also plotted is the sample of observed disks that were 
actually detected at 24 and 70 ^m (§4.2.1). To a large 
degree the observed population is well described by the 
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Fig. 4. — Fractional luminosity versus radius. (Left) The 10,000 disks in the A star population model of Fig. 2 are shown with dots. 
The detection thresholds of the Su et al. (2006) Spitzer observations at 24 and 70 lira (i 7, 24disk/^ 7 24* > 0.1 and -FVodisk/^TO* > 0.55) arc 
shown with solid lines for A0V stars (equation 19). The maximum luminosities for planctcsimal belts around A0V stars of age 8, 80, 280, 
and 800 Myr in the model are shown with dashed lines (equation 20). (Right) As (left), except that only disks in the model population 
that could be detected at 24 and 70 fim are plotted. Also shown with filled circles are the properties of the 46 disks that were detected at 
both 24 and 70 (im (§4.2.1). The dotted line shows the limit below which P-R drag becomes important for A0V stars (equation 21). 



model population, both in terms of the lower envelope 
(described by the 24 and 70 /zm detection limits), and 
the upper envelope (described by a oc r 7 / 3 power law and 
close to the line of maximum luminosity for disks a few 
Myr old), as well as the distribution of disks in between, 
which is more dense between the lines of maximum lu- 
minosity for 80 and 800 Myr. A fit to the population of 
observed disks gives / oc r -° 1±0 - 3 with a large spread 
at each radius, which is not inconsistent with that ex- 
pected from the model population for disks detected at 
24 and 70 /zm (/ oc r 6 ). The agreement is closer (fit to 



observed population of / 



„0.9±0.3 



) when the systems 



mentioned in the following paragraph are removed. 

Closer inspection shows that the disks that lie above 
the 8 Myr maximum luminosity line for disks around 
A0V stars are: the disks with well known ages in 
the range 5-12 Myr (HD39060, HD109573, HD1 10058, 
HD141569 and HD172555), and three stars that have 
///max in the range 30-100 and so appear unusually 
bright for their age (HD3003 at 50 Myr, HD38678 at 230 
Myr, and HD115892 at 350 Myr; Table 1). The pres- 
ence of the disks of 5-12 Myr stars above the 8 Myr line 
is to be expected, since the planetesimal belts in these 
systems either must be, or are likely to be < 8 Myr-old. 
Even if the star itself is older than 8 Myr, we still need 
to allow some time for the planetesimal belt to form, 
and if this occurs on the same timescale that it takes 
for proto-planetary disks to dissipate, then it takes ~ 6 
Myr (Haisch, Lada, & Lada 2001). Also, at such a young 
age, there is the possibility that these systems have yet 
to reach equilibrium, and that this early phase is char- 
acterized by an unusually high quantity of small dust 
grains. Whether the three unusually bright stars are in 
a transient phase (Wyatt et al. 2007), or have unusual 
planetesimal belt properties remains to be seen. 

Understanding how the detection limits (equation 19) 
and constraints from the maximum possible luminosity 
for a given age (equation 20) appear on the plot of frac- 
tional luminosity against radius (left Fig. 4) readily ex- 
plains trends (i) to (iii) given in §4.2 (Su et al. 2006): (i) 
The model predicts that disks that are detected at 70 /im 
but not 24 /zm should lie in the region below the 24 /zm 



threshold and above the 70 /zm threshold on this plot, 
with / in the range 10~ 6 to 3 x 10~ 5 (for A0V stars), 
depending on their radius. In contrast, those detected at 
both wavelengths, which should lie above both thresh- 
olds, would have / in the range 5 x 10~ 6 to 10~ 2 (for 
A0V stars). Thus the model predicts that the range in 
/ should be systematically higher for disks detected at 
both wavelengths than for those detected only at 70 ^m, 
as found by Su et al. (2006). (ii) The region of param- 
eter space occupied by disks that are detected at 70 /zm 
but not at 24 /zm also lies entirely below the line of max- 
imum luminosity for 280 Myr (for A0V stars; Fig. 4). 
This means that the only disks in this region of parame- 
ter space that are younger than 280 Myr must be those 
of initially low mass for which the most massive objects 
within them have yet to come to collisional equilibrium. 
This means that this population has a higher average age 
(460 Myr) than the rest of the population (400 Myr). 
This is consistent with the findings of Su et al., since the 
mean age of their group II sources (i.e., those detected at 
70 but not 24 /zm) is higher than average at ~ 490 Myr. 
The model can also be used to analyze the population of 
disks that are expected to be detected at 24 but not 70 
/zm (i.e., Su et al. 2006 group V sources). In the model 
such disks have small radii (< 23 AU for A0V stars) and 
lie above the 280 Myr line of maximum luminosity (for 
A0V stars), meaning that there can be no disks in this 
population that are older than this limit, and that the 
mean age of this population is ~ 150 Myr. Indeed, the 
oldest star in Su et al. (2006) group V is 400 Myr old, and 
the mean age of their group V sample is 100 Myr. 5 (iii) 
The lowest excess in the model population is / w 10~ 9 , 
and 28% of stars in this population have excesses lower 
than 10~ 7 , in agreement with non-detections at this level 
in Su et al. (2006). 

4.2.3. Fractional luminosity vs age 

5 Some of the Su et al. group V sources were not observed down 
to the limit -f?70dct = 0.55, e.g., due to cirrus confusion. Thus more 
sensitive observations may find 70 (im excesses at i?70det > 0-55 
in this sample. However, this would not affect the conclusion that 
bona-fide group V sources are, on average, younger than the rest 
of the population. 
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Fig. 5. — Fractional luminosity versus age. (Left) The 10,000 disks in the A star population model of Fig. 2 arc shown with dots. 
The evolutionary tracks of individual planetesimal belts at 10 AU and 100 AU around A0V stars are shown with solid and dashed lines, 
respectively, for starting disk masses of 1, 10 and lOOMgj (with higher masses corresponding to higher fractional luminosities at young 
ages). (Right) As (left), except that only disks in the model population that could be detected at 24 and 70 fim are plotted. Also shown 
with filled circles are the properties of the 46 disks that were detected at both 24 and 70 (im (§4.2.1). The solid and dashed lines show fits 
to the observed and model disk populations, respectively. 



The fractional luminosity of the 10,000 disks in the 
model population are shown in left Fig. 5 as a function 
of stellar age. Also shown on this plot are lines show- 
ing the evolution that a planetesimal belt at 10 and 100 
AU would take assuming starting masses of 1, 10 and 
100M®. As described in §2.1, this evolution remains flat 
until the largest planetesimals collide at which point the 
fractional luminosity falls off oc t~^ e . For the planetesimal 
belt at 100 AU, the largest planetesimals do not have a 
chance to reach collisional equilibrium over the assumed 
800 Myr main sequence lifetime of the star, while for the 
planetesimal belt at 10 AU, this equilibrium is reached 
in the range a few Myr to a few hundred Myr, depend- 
ing on starting mass. The upper envelope in the model 
population can be described as being reasonably flat at 
/ w 10~ 3 until an age of ~ 100 Myr at which point it falls 
off. The fall-off rate is slightly flatter than oc t~^ e , which 
is because the brightest disks at any given age are at dif- 
ferent radii for different ages. For example, the upper 
envelope for ages < 100 Myr is reasonably well charac- 
terized by the evolution of the lOOM® planetesimal belt 
at 10 AU shown on left Fig. 5, whereas the evolution 
of the IOOMq planetesimal belt at 100 AU has a higher 
luminosity than that at 10 AU at ages > 600 Myr. 

The sub-sample of the model population that can be 
detected at 24 and 70 ^m is shown in right Fig. 5 and 
exhibits the same upper envelope. The lower envelope of 
this population also falls off slowly with time, not because 
the detection limit drops, but because a smaller fraction 
of the stars at younger ages have luminosities approach- 
ing that limit, whereas collisional processing has reduced 
the fractional luminosities of older stars so that many 
are pushed below the detection threshold (particularly 
for disks at small radii). The model thus explains trends 
(iv) and (v) of §4.2 (Su ct al. 2006): (iv) The spread 
in fractional luminosity at any given age is around two 
orders of magnitude both in the observed and model pop- 
ulations, (v) Both the model and observed populations 
have lower fractional luminosities for disks around older 
stars than around younger ones, with the model popu- 
lation exhibiting a / oc i age 39 fall-off that is in excel- 
lent agreement with the observed population for which 



/ *ago' 57±0 ' 13 (compared with / oc i ag c' 6 quoted in Su 
et al. 2006 for their observations). 

The biggest difference between the model and observed 
populations is that the model assumes that stellar ages 
are uniformly distributed between and 800 Myr, which 
means that the model predicts that 30% of the popula- 
tion detected at 24 and 70 /im should be older than 400 
Myr, whereas just 11% of the observed population are 
that old. The most likely cause for this discrepancy is 
the fact that the mean age of the stars searched for 24 and 
70 /im excesses is younger than the 400 Myr mean age 
of the model population. The mean age of the Su et al. 
(2006) sample is 270 Myr which may be attributable to a 
bias in the way the sample was chosen, or to the fact that 
A0V stars reach the end of the main sequence at ~ 400 
Myr which means that a dearth of > 400 Myr early type 
stars is to be expected in any sample of main sequence 
stars (Greaves & Wyatt 2003). The fact that the model 
provides a close fit to the fraction of stars from any given 
age range that lie in specific ranges of fractional excess 
flux (§3.2 and §4.1) means that the age distribution of 
the observed sample is immaterial suggesting that a sud- 
den disappearance of disks at ~ 400 Myr (Habing et al. 
1999) is not required to explain the statistics. However, 
given that 39 stars in Su et al. (2006) were observed in 
the > 400 Myr population, and that the model predicts 
that 46% (i.e., 18) of these should have disks above the 
nominal 70 /am detection threshold quoted in that pa- 
per, then the observed number of detected disks of 3 is 
only consistent with the model when one takes into ac- 
count that not all of the 39 stars were observed down to 
-fiVudet = 0.55, e.g., because of cirrus confusion. Almost 
all (35) of the 39 were observed to i?7odot = 5, and the 
model predicts that 17% (i.e., 6) of these should have 
been detected at this level, which is consistent with the 
observed number. In other words we predict that, unless 
there is some mechanism destroying these disks at ~ 400 
Myr, then more disks remain to be detected around the 
> 400 Myr star sample at the nominal 70 /zm detection 
threshold. 

The evolution of the upper envelope discussed in trend 
(vi) (§4.2; Su et al. 2006) is also in agreement with the 
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model, since we find a fall-off in the upper envelope that 
is close to oc t~g C 9 . In answer to the discrepancy with 
the Decin et al. (2003) result that showed a maximum 
fractional luminosity that is constant with age, we first 
note that this applied to disks around stars of all spectral 
types. When considering their results only for A stars, 
we find an evolution which is in agreement with that of 
our observed sample, with the exception of HD22128, an 
A5 star for which Decin et al. (2003) quote an age of 
<~ 1.4 Gyr and a fractional luminosity of 70 x 10~ 5 (i.e., 
which would lie off to the right of our right Fig. 5 some- 
what above the predicted model population). We thus 
consider that the anomalous system HD22128 could ei- 
ther have disk properties that make it unusually bright, 
it could be in a transient state, or its luminosity or age 
could have been miscalculated. Indeed an age of 320 
Myr has also been quoted for HD22128 (Decin et al. 
2003). Further, recent MIPS observations show that the 
source is extended at 24 and 70 /im on > 1000 AU scales. 
Since such emission is unlikely to originate in a disk, this 
suggests that the luminosity of the debris disk based on 
bolometric far-IR observations is an upper limit and the 
existence of a debris disk in this system is suspect. 

4.2.4. Radius vs age 

The radii of the 10,000 disks in the model population 
are shown in Fig. 6 as a function of stellar age. Since in 
this model the radius does not change with age we have 
plotted only those disks in the model that it is possible 
to detect at 24 and 70 /j,m and compared these with the 
properties of the observed population. The most notice- 
able trend of the model population is that at any given 
age it is only disks that have radii larger than a certain 
limit that can be detected, and that this limit increases 
with age. This can be readily understood from left Fig. 4. 
For example, disks that are both <~ 80 Myr old and de- 
tected at 24 and 70 /xm must lie below the 80 Myr line 
of maximum luminosity on that figure, and above the 
24 and 70 fim thresholds. This is a region of parame- 
ter space that permits only disks with radii > 16 AU 
(for A0V stars), and that does not extend down so far 
in radius at later ages (and encompasses a much smaller 
area). Physically this is because planetesimal belts at 
smaller radii process their mass faster and so fall below 
the detection threshold at younger ages (i.e., it does not 
mean that planetesimal belts at this radius do not ex- 
ist). This is manifested as an apparent increase in the 
radius of planetesimal belts with age, with a least squares 
fit to the population of detected disks in the model in- 
creasing oc i° g o°, although there is still a large spread in 
radii at any age. While the effect is small, the increasing 
size of the empty region with age means that the radius 
distribution (and so the color temperature distribution) 
is broader at younger ages, which means that the model 
also reproduces trend (vii) (§4.2) seen by Su et al. (2006). 
The modest apparent increase in size with age is also con- 
sistent with the observed population, for which we derive 
a reasonably flat evolution of oc i" g ° 6±0 - 07 (although note 
that a steeper evolution oc iago 9±0 '° 5 would be inferred 
if the anomalous sources discussed below are removed). 
A flat distribution of radii with age was also observed 
in the sub-mm study of Najita & Williams (2005) which 
included disks around stars of all spectral types. It was 
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Fig. 6. — Planetesimal belt radius versus age. The model pop- 
ulation of Fig. 2, shown with dots, includes just those disks that 
could be detected at both 24 fim or 70 (im. Also shown with filled 
circles are the properties of the 46 disks that were detected at 
both 24 and 70 fim (§4.2.1). Least squares fits to this evolution 
are shown with solid and dashed lines for the observed and model 
populations, respectively. 

also reported by Rhee et al. (2007), although these au- 
thors suggested that the upper envelope in radius may 
increase across the range 100 to 1000 Myr, an observa- 
tion which is not predicted in our model, and only seen 
with low significance in our observed population. 

However, there is a distinct anomaly in that six of 
the observed disks fall in the region of Fig. 6 in which 
there are no model disks: HD38678, HD1 15892, and 
HD172555 lie well within the forbidden region, while 
HD2262, HD3003, and HD 106591 lie close to the upper 
edge of this region. We offer two possible explanations 
for these systems. The first is that these are systems 
which are undergoing a transient period of high dust lu- 
minosity (Wyatt et al. 2007). This has already been 
discussed as a possibility for HD38678 in §4.2.2 based on 
its high luminosity for its age (and was also suggested in 
Moerchen et al. 2007). The other possibility is that the 
fitting procedure has underestimated the radius of the 
planetesimal belt. 

In addition to the possibility that a dust temperature 
that is higher than black body has caused the dust loca- 
tion to be underestimated, it is also possible that the dust 
is located closer to the star than the planetesimal belt be- 
cause of the action of Poynting- Robertson (P-R) drag. P- 
R drag becomes important when 10 4 } '(r / 'dr) yJr/Ml / ' f3 < 
1 (Wyatt 2005), where [3 is the ratio of the radiation force 
to that of gravity of the star, since this is when the col- 
lisional lifetime of the dust is equal to the time it would 
take to reach the star by the drag force. Assuming, as in 
the model here, that the dust luminosity is dominated by 
the smallest dust in the distribution with [3 — 0.5, this 
limit is given by / < / pr , where 



/ pr - 50 x 10- 6 (dr/r)^/M+/r, 



(21) 



and this is shown on right Fig. 4 for A0V stars assuming 
dr/r = 0.5, and compared with the observed fractional 
luminosity for the 46 disks in Table 1. The proximity of 
this P-R drag limit to the detection threshold at 24 and 
70 /jm illustrates how P-R drag may become important 
for disks close to the detection threshold of Spitzer, at 
least for early type stars. 
The details of how the model would be affected by P-R 
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Fig. 7. — Histogram of the infrared luminosity divided by the 
maximum luminosity possible in the model for a disk of this age, 
spectral type and radius (assuming planetesimal belt properties 
inferred in §4). The distribution for the observed sample of 46 disks 
that were detected at both 24 and 70 fim (§4.2.1) is shown with 
a solid line. The distribution for disks in the model population of 
Fig. 2 that could be detected at 24 and 70 fim scaled to a population 
of 46 disks is shown with a dashed line. A Gaussian smoothing of 
1 dex applied to the model population distribution is shown with 
a dotted line. 

drag below (or approaching) this level are not clear with- 
out further modeling, but the nature of the drag force 
means that a lower dust radius (but not planetesimal belt 
radius) would be expected. Thus this is one possible ex- 
planation for the close proximity of dust from the stars 
HD2262, HD 106591, and HD 11 5892, all of which have 
///pr < 1. It may be possible to test this observation- 
ally, since in the extreme situation where P-R drag dom- 
inates the dust distribution, the surface density of the 
dust disk would be uniform from the planetesimal belt 
all the way in to the star leading to an emission spec- 
trum that increases linearly with wavelength (F u oc A; 
Wyatt 2005). One other source for which P-R drag may 
be important in shaping the inner edge of its dust disk is 
HD 19356, for which /// pr < 1, while a further 11 disks 
have 1 < /// pr < 5 (Table 1). 

4.2.5. Accuracy of / max 

Figure 7 shows the ratio of observed fractional infrared 
luminosity to the maximum luminosity expected for disks 
at the inferred radius around stars of the age and spec- 
tral type of the hosts; i.e., ///max- This is found for the 
sample of 46 stars detected at 24 and 70 fim (§4.2.1) to 
be peaked at ///max ~ 1, but extending to around a 
factor of 100 higher and lower; the values for this sam- 
ple are also given in Table 1. The sub-sample of the 
model population detectable at 24 and 70 ^tm is also 
shown on this figure, and this peaks much more sharply 
at ///max ~ 1, with a tail to lower ///max- As expected 
there are no disks in the model with / > / max , and the 
fact that the majority of the model disks have / w / max 
indicates that the majority of the detected disks (in the 
model) are evolving with their largest planetesimals in 
collisional equilibrium. 

The observed population has a noticeably broader dis- 
tribution of /// ma x than does the model population. We 
interpret this to be a consequence of the assumption that 
all disks have the same planetesimal properties. In re- 
ality disks would be expected to have a distribution of 
properties and the parameters (D c , and e) used in the 



expression for / ma x would only be representative of the 
population as a whole rather than of individual disks. 
Thus the breadth of the observed distribution may be 
indicative of the distribution of disk properties. To esti- 
mate this breadth we show on the same figure a smoothed 
model in which the distribution has been smoothed by a 
Gaussian of width 1 dex, which fits the observed distri- 
bution well, and suggests a similar breadth in the distri- 
bution of Z?°- 5 <5p 5 / 6 e -5 / 3 about the mean values given 
in equation (17). Since the ///max distribution will also 
have been smeared out due to uncertainties in the ob- 
served (or inferred) disk properties, notably r and /, the 
true distribution of £>°- 5 <2p 5 / 6 e~ 5 / 3 is likely to be much 
smaller than 1 dex. It is, of course, also possible that 
some of the disks inferred to have / > /max are m a 
transient phase (Wyatt et al. 2007), which would sug- 
gest an even narrower distribution of D®- 5 Q^ 5 / 6 e~ 5 / 3 for 
the population of disks that are in steady-state. This in- 
dicates a large degree of uniformity in the properties of 
A star disks. 

The most likely candidates for transient emission 
around A stars are readily identifed in Table 1 as those 
with ///max > 10: HD3003, HD38678, HD1 15892, and 
HD172555. However, since none of these reaches the 
threshold /// max > 100, we do not claim that any of 
these must be transient, but point out that their prop- 
erties would have to be extreme compared with the rest 
of the A star population for their dust emission to be 
explained by steady-state collisional processing of plan- 
etesimal belts. 6 The inferred breadth of the distribution 
indicates that a disk with /// max > 100 is a 2cr devi- 
ation, while ///max > 1000 is a 3a deviation. We also 
note that the black body assumption would tend to over- 
estimate the value of ///max, since underestimating the 
radius by a factor of 3 would lead to an overestimation of 
/max by a factor of 13 and so reduce the factor /// max by 
an order of magnitude. In other words, it is not possible 
to conclude that the disks with /// max = 10— 100 repre- 
sent a significant deviation from the disks found around 
the remainder of the A star population. 

5. PREDICTIONS FOR SCUBA-2 LEGACY SURVEY 

While the fact that the model population can explain 
the available far-IR statistics with steady-state collisional 
evolution is important, this model population can also 
be used to make predictions for what we can expect to 
find in future debris disk surveys. Since the model was 
originally constrained by a fit to the 24 /im statistics, 
in essence the first of those predictions was for the far- 
IR surveys at wavelengths other than 24 /xm, which §4 
showed to be confirmed. This section demonstrates the 
further applicability of the model, by making predictions 
for the outcome of the SCUBA-2 legacy debris disk sur- 
vey. This is an unbiased 850 /im survey of the nearest 
100 stars in each of the spectral types A-M (i.e., 500 stars 
in total) (Matthews et al. 2007), and will be sensitivity 
limited down to the confusion limit which gives a 3cr limit 
of ^ssoiim = 0.002 Jy. For A stars the survey extends out 
to around 45 pc. 

6 Note that a size of 3 AU for HD38678 as inferred from mid- 
IR imaging (Mocrchcn et al. 2007) would reduce its maximum 
possible fractional luminosity given its age, /max, by a factor of 7, 
and so would support the transient nature of this source. 
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The model is able to make predictions for the 850 /tm 
flux of the A star population through application of equa- 
tion (18) to get the disk flux, as well as using equation (8) 
to get the stellar flux (although this generally falls below 
the detection threshold). Comparing the flux predicted 
by equation (18) with that of the 7 A stars in the litera- 
ture with measured 850 /im fluxes (HD109573, Greaves, 
Mannings & Holland 2000; HD39060, HD216956, and 
HD172167, Holland et al. 1998; HD141569, Sheret et 
al. 2004; HD21997 and HD14055, Williams & Andrews 
2006) shows that this equation overestimates the disk 
flux by an amount X^o, where X^5o for this population 
lies in the range 2-8, with a median value of 4.0. This is 
to be expected, since it is known that the emission spec- 
trum falls off steeper than a black body at wavelengths 
beyond ~ 70 /im (e.g., Dent et al. 2000), because of the 
low emission efficiency of the small dust grains which 
dominate the cross-sectional area in the disk. The factor 
X 8 5o is expected to depend on a combination of the dust 
size distribution, compositional properties, planetesimal 
belt radius, and stellar spectral type. Thus, when deriv- 
ing 850 /im fluxes for the model population we reduced 
the flux found from equation (18) by a factor of Xg5o = 4. 
The ability to estimate sub-mm flux with reasonable ac- 
curacy from extrapolation of far-IR data is in agreement 
with the conclusion of Rhce ct al. (2007) who found a 
strong correlation of dust masses derived from sub-mm 
fluxes, M S50 = F v d 2 [n l/ B l/ (X, T)] _1 , with / derived from 
far-IR fluxes. Indeed, our model reproduces the trend 
shown in their Fig. 5, since with a dust opacity of 0.17 
m 2 kg~ x and a factor of A 850 = 4 our model predicts 
//M 850 = 14r~ 2 (in units of M^ 1 ), exactly as observed 
by Rhee et al. (2007) for A stars. 

Figure 8 shows the predictions of the model for the 
fractional luminosity versus radius of the sub-sample of 
the model population that have 850 /im fluxes above the 
detection threshold of 2 mJy (noting that all of these are 
also within the survey threshold of 45 pc). Also shown 
are the 7 A star disks with sub-mm detections, indicating 
the 4 of these that are within 45 pc. All but one of 
the 4 within 45 pc will appear in the SCUBA-2 survey; 
HD39060 is excluded, since it is below the declination 
limit of —40°. There is good agreement between the 
model and observed populations inside 45 pc, insofar as 
it is possible to say with such small number statistics. 
Since the sub-mm detection limit is determined by the 
instrumental sensitivity, rather than by the accuracy of 
the photospheric calibration, it is no longer described by 
equation (19), but can be given by 

/ dot = 4.3xl0 10 A 8 5o J F 1 850iim(rf/r) 2 /^(A,278.3i^ 25 r- - 5 ), 

and so is distance dependent. This limit falls with plan- 
etesimal belt radius, because for disks at a given distance 
it is easier to detect those of lower fractional luminosity if 
they are at larger distance from the star, since their lower 
temperature at larger distances is more than compen- 
sated by the larger cross-sectional area of material. The 
SCUBA-2 detection limit is plotted on Fig. 8 for A0V 
stars at 5 and 45 pc assuming their disks have A 850 = 4. 
The lower envelope of the model population is thus nat- 
urally explained by the small number of A stars that are 
close to us (of order 14 of the 10,000 are within 5 pc). It 
is tempting to suggest that the lower envelope of the dc- 
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Fig. 8. — Predictions of the model population of Fig. 2 for the 
outcome of the SCUBA-2 legacy survey, an unbiased survey at 850 
(im of the nearest 100 A stars. (Top) Fractional luminosity of 
the sub-sample of the model population which is both within the 
survey limits (< 42 pc), and which has a predicted flux above the 
sensitivity limit of the SCUBA-2 survey (-Fgsodisk > 2 mJy) (dots). 
The solid lines are the detection thresholds of 24 and 70 fim Spitzcr 
observations of Su et al. (2006), and the 850 fim detection thresh- 
old for disks at 5 and 42 pc in the survey; all of these detection 
thresholds apply to disks around A0V stars. The disks of A stars 
at < 42 pc that have already been detected at 850 fim are shown 
with filled circles, while those that have been detected but lie be- 
yond 42 pc are shown with diamonds. (Bottom) Histogram of the 
distribution of 850 fim flux expected for the 100 A stars within 42 
pc. The detection limit is shown with a dotted line. The dashed 
line shows the stars in the survey already detected at 850 fim. 

tcctcd disks is explained in the same way, although with 
just 4 disks the lower envelope is not so meaningful. In 
any case, it is the nearest disks which are detectable in 
the sub-mm, while for those in a given distance range it 
is predominantly large disks, or those closer in with high 
fractional luminosity, that can be detected. 

Bottom Fig. 8 shows the histogram of the 850 /tm 
fluxes of the model population (noting that this popu- 
lation should be representative of the SCUBA-2 survey 
sources, since all are within 45 pc). The model predicts 
that 17 of the 100 stars in the SCUBA-2 A star sam- 
ple should have a flux above the detection threshold of 
2 mJy. Of those 17: five would be expected above 10 
mJy, and eight in the range 3-10 mJy, with the remain- 
ing four detected with 3 — 4.5cr confidence at 2-3 mJy. 
The histogram also shows the sub-mm flux distribution 
for the 3 A stars already detected, showing that we al- 
ready know of 2 stars within 45 pc with fluxes above 10 
mJy (HD172167, HD216956), while we know of just one 
in the range 3-10 mJy (HD14055). Thus the SCUBA-2 
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survey is expected to discover three new bright disks > 10 
mjy, but can expect to find a further seven new disks at 
3-10 mjy levels. The > 10 mjy sources are particularly 
interesting, since at such a level of flux it is possible to 
resolve the structure of the disks (Holland et al. 1998; 
Wyatt et al. 2005), and indeed the model predicts that 
60% of these will have an angular size of > 6 "(noting 
that all may be larger than this given that the black body 
assumption underestimates the true disk size by a factor 
of up to three, Schneider et al. 2006). It is also possible 
to extract some information on the spatial distribution 
of the 3-10mJy sources, e.g., giving a direct measure of 
the size and inclination of the disk (Greaves et al. 2004) , 
with around 23% of these predicted to be larger than 6 

It may even be possible to extract information from 
the 83 stars without detected disks, since they are pre- 
dicted to have a mean flux of 0.60 mjy, of which 0.26 mjy 
is the mean disk flux and 0.34 mjy is the mean stellar 
flux. By co-adding their data, the statistical noise level 
would be reduced to 0.7/-\/83 = 0.077 mjy, suggesting 
that it should be possible to get a positive detection of 
the photospheric emission, and also to discern if there is 
an excess coming from the disks with <~ 3a confidence. 
However, contamination from extragalactic sources that 
lie within the SCUBA-2 beam would have a signal of 
comparable magnitude to the mean disk flux and would 
reduce our ability to make a conclusive detection of disk 
emission from this population. 

The model predicts that there will be little age depen- 
dence within the sample of disks detected by SCUBA-2, 
in agreement with the result of Rhcc et al. (2007) who 
estimated the mass evolution based on far-IR observa- 
tions. For example, the model predicts that the mean 
age of the detected sample is ~ 340 Myr (compared with 
400 Myr for the whole population), with a roughly even 
number of disks (8-9) detected in both the < 300 Myr 
and > 300 Myr age bins, which will be indistinguishable 
from that expected if there were no age dependence (6 
in the < 300 Myr and 11 in the > 300 Myr samples). 
Naturally SCUBA-2 is inclined to detect disks with low 
values of Xg,^, i.e., those for which their emission is more 
like black body (which would be inferred to be composed 
of large grains). However, given the small range in A§5o 
observed, this dependence is very slight and would af- 
fect mainly those disks detected close to the sensitivity 
threshold. There is also predicted to be a small spectral 
type dependence, in that 19% of stars later than A3V 
are predicted to be detected by SCUBA-2, compared 
with 15% of those earlier than A3V. This arises because, 
while disks of lower fractional luminosity can be detected 
around earlier type stars, /dot L~ ' 25 (equation 22), 
the fractional luminosities of their dust belts are also 
much lower because of the higher dust blow-out radius, 
/max oc L~°- 5 (equation 14). We do not anticipate that 
this would be detectable within the A star sample, how- 
ever, it suggests that the detectability of sun-like stars 
could be high, should their disk population be described 
by the same parameters as for the A stars. 

Top figure 8 also shows the 24 and 70 /im detection 
limits of the Spitzer survey of Su et al. (2006). This 
shows that almost without exception (99.9%) the disks 
that can be detected by SCUBA-2 at 850 ^im are also 



detectable by Spitzer observations at 70 fim with a de- 
tection threshold of i?7odot > 0.55, while most (61%) 
of those detectable by SCUBA-2 are also detectable by 
Spitzer at 24 /im. In other words, ten of the 100 A stars in 
the SCUBA-2 sample should be detectable at 24, 70 and 
850 /im, while seven should be detectable at 70 and 850 
/im, but not at 24 /im; the combination of sub-mm and 
far-IR data would be very important for constraining the 
dust temperature in these disks. While this suggests that 
the survey could equally well be done at 70 /<m, there are 
three important points to consider, (i) First, not all of 
the sample has been observed with Spitzer - the Su et 
al. (2006) sample includes 42 sources within 45 pc (of 
which 20 were detected at 70 /im), which assuming even 
coverage across the sky means that 34 of these should 
fall in the SCUBA-2 survey (which has declination lim- 
its of -40° and +80°). This implies that two thirds of 
the SCUBA-2 A star sources have not previously been 
searched for disks in the far-IR (although bright disks 
may have turned up in IRAS surveys which have a de- 
tection threshold well above that shown in Fig. 8). It is 
also not possible to scrutinize all stars with Spitzer to 
the level of i?7odct = 0.55 due to cirrus confusion. Thus 
it is predicted that the SCUBA-2 survey will turn up a 
significant number of genuinely new discoveries, (ii) Sec- 
ond, it is important to emphasize the unbiased nature of 
the SCUBA-2 survey. Several sources have recently been 
identified as possessing excess emission in the sub-mm 
even when no excess is present in far-IR surveys (Wy- 
att et al. 2003; Holmes et al. 2003; Najita & Williams 
2005). These disks must be cool (< 40K) and at large 
radius; e.g., extrapolation of the detection thresholds in 
Fig. 8 to larger radii indicates that the sub-mm is more 
sensitive than 70 /im surveys to disks that arc > 210 AU 
at 5pc for A0V stars and > 74 AU at 5 pc for A9V stars. 
Since the model population was constrained by 24 and 
70 /im surveys, there was no need to invoke a popula- 
tion of disks of large radius, although we know that dust 
exists around young A stars to hundreds of AU (Kalas 
& Jewitt 1995; Clampin et al. 2003). Thus the pre- 
dicted detection statistics should be viewed as a lower 
limit, with many large disks being detected at 850 /im 
which cannot be detected in the far-IR. Previous unbi- 
ased sub-mm surveys suggest that ~ 15% of stars could 
have disks too cold to detect in the far-IR, which could 
push the detection rate to - 30% for the SCUBA-2 A 
star survey, (iii) Even discounting the ultra-cool disks 
discussed in (ii), the sub-mm statistics will set important 
constraints on the distribution of the planetesimal belt 
radii, since the predictions made here are very sensitive 
to that distribution. For example, with a radius distribu- 
tion N(r) cx r -0 ' 5 , which is not ruled out in the present 
study at the la confidence level, the model would predict 
that 22 of the 100 A stars would be detected at > 2 mjy 
of which 7 would be > 10 mjy. 

6. CONCLUSIONS 

A simple model for the steady-state evolution of dust 
luminosity for planetesimal belts evolving due to colli- 
sions was described in §2. This section also described 
how the model could be applied to determine the prop- 
erties of the disks of a population of stars, given that 
these disks would have had a range of initial properties. 
This was applied in §3 to the population of A stars that 
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was searched for evidence of dust emission at 24 /zm by 
Rieke et al. (2005), and it was shown that their detec- 
tion statistics as a function of age could be reproduced 
with a model population in which the largest planetes- 
imals have a size of D c = 60 km, and all planetesimals 
have a strength of QJ, = 300 J kg -1 and an eccentricity 
of e = 0.05. However, these parameters should not be 
over-interpreted, since the observable properties of the 
model population remain the same as long as the com- 
bination Z?"- 5 <5p 5 ^ 6 e -5 / 3 is unchanged. Thus more de- 
tailed models of the collisional evolution of planetesimal 
belts (e.g., Krivov, Sremcevic & Spahn 2005) are needed 
to interpret these parameters, e.g., in terms of the com- 
position of the planetesimal belts. For now we note that 
these are reasonable parameters based on planetesimal 
(Benz & Asphaug 1999) and planet formation models 
(Kenyon & Bromley 2002). The model used as input 
a distribution of planetesimal belt starting masses that 
is log-normal with the same width (la of 1.14 dex) as 
inferred for proto-planetary disks (Andrews & Williams 
2005), and found such a distribution would have to be 
centred on M m id = lOM®. This mass is consistent with 
that expected for proto-planetary disks around A stars, 
although the observable properties of the model popu- 
lation are also reproduced as long as M m iaD~ 5 is un- 
changed. In other words, the observable 24 /j,m prop- 
erties of main sequence A star disks can be explained 
by the steady-state evolution of a disk population with 
realistic starting parameters (although more extreme so- 
lutions are also possible). 

In §4 this model was tested against its predictions for 
the 70 /zm properties of the disks in this population. 
These were found to reproduce the statistics of the Su 
et al. (2006) survey which showed a much longer de- 
cay timescale for the 70 fim excess than for that at 24 
(Um. The model was also used to predict the properties 
of the sub-sample of disks that could be detected at both 
24 and 70 /zm, and this was compared with a sample of 
46 A stars compiled from the literature for which excess 
emission has been detected at both 24 and 70 fim. The 
model reproduces and explains the distribution of the 
observed disks on the / vs r (Fig. 4), / vs t agc (Fig. 5) 
and r vs i age (Fig. 6) plots, including an upper envelope 
in the f vs r plot that increases oc r 7 ' 3 , a decay in the 
mean luminosity of disks with age oc i a g e ' , and a mean 
radius of detected disks that, if anything, increases with 
age oc t° g ^°. 

Thus it appears that the wide range in typical 
disk properties can be explained without appealing to 
stochasticity. The large spread in 24 /zm fractional ex- 
cess at any age (Rieke et al. 2005) can occur naturally 
from a spread in initial disk masses and radii, and the 
slower fall-off at 70 /zm and more detailed distributions 
of the properties of disks detected at 24 and 70 ztm (Su et 
al. 2006) could be the consequence of the detection bias. 
This bias is best explained on the plot of / vs r, on which 
it is possible to plot the lines above which disks must lie 
to be detected at 24 and 70 /zm, and below which they 
must lie for a given age. One illustrative example of this 
bias is the fact that while in the model the planetesimal 
belt of any individual star does not change in radius, the 
mean radius of detected planetesimal belts is predicted to 
increase with age. The reason is that planetesimal belts 



that are closest to the star are processed faster and so fall 
below the detection threshold much faster than planetes- 
imal belts that are farther from the star. Thus if such 
a trend is seen it would not necessarily be evidence for 
delayed stirring, in which the collisional cascade is only 
initiated once enough time has elapsed for Pluto-sized ob- 
jects to form, a process which takes longer further from 
the star (e.g., Kenyon & Bromley 2002; Dominik & Decin 
2003). Rather it could be evidence that the processing of 
the planetesimal belt is occurring inside-out, as inferred 
by Su et al. (2006). 

This does not rule out the possibility that some, or even 
all, A star debris disks in the Rieke et al. (2005) and Su 
et al. (2006) samples are undergoing transient events on 
top of the overall trend established by the model. In 
fact, there are a few disks that may (rather than must) 
require this mechanism to explain their properties. The 
potentially transient systems identified in this paper are: 
HD3003, HD38678, and HD172555, which have high lu- 
minosity for their age; as well as HD1 15892, which has 
unusually an low radius for its age. All of these sys- 
tems are identifiable by their high value of /// max , where 
/max is the maximum possible fractional luminosity that 
a disk can have given its age and radius (assuming the 
planetesimal belt properties inferred for the rest of the 
population). We also found that the disks of HD2262 
and HD106591 have unusually low radii for their age, 
and suggested that this may be attributable to the ac- 
tion of P-R drag on dust produced in their dust belts, 
since this starts to become important when / < / pr , and 
these belts (along with HD115892 and HD19356) can be 
identified from the sample as having /// pr < 1- Neither 
do the results in this paper imply that the collisional 
cascade in all A star planetesimal belts must have been 
initiated soon (~ 10 Myr) after the stars formed, rather 
than the collisional cascade being initiated after a delay 
of several 10s or 100s of Myr in some systems (i.e., the 
delayed stirring model) . Indeed, if the lack of large radii 
disks at young ages is confirmed (Rhee et al. 2007), this 
would support some role for delayed stirring in debris 
disk evolution. 

Analysis of the distribution of ///max of the 46 A star 
disks detected at both 24 and 70 fim shows a distribu- 
tion that may be attributed to the planetesimal belts 
not having exactly the same properties, but having a 
distribution in the combination of parameters given by 

Dc' 5 Qr> 5 ^ e ~ 5 ^ 3 (which is also the combination in the 
parameter /max) that is log-normal with a la width of 1 
dex. This means two things: First, that the distribution 
of parameters of these A stars is quite uniform. This is 
perhaps to be expected if the planetesimals in all belts 
are of similar composition, and grow to a similar maxi- 
mum size before perturbations from the largest member 
of the cascade (around 2000 km) stir the rest of the plan- 
etesimal belt resulting in eccentricities of similar magni- 
tude (e.g., Kenyon & Bromley 2002). Second, that the 
anomalous systems identified in the last paragraph are 
only l — 2cr anomalies in that they may be explained as 
disks with higher than average planetesimal strength or 
size, or lower than average eccentricity. In contrast, this 
analysis strengthens the result of Wyatt et al. (2007), in 
which the dust emission from several sun-like stars, such 
as HD69830, was inferred to be transient, since if the 
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planctesimal belts of sun-like stars are described by sim- 
ilar properties to those of A stars, then all of the disks 
they inferred to be transient have ///max > 4000 and 
each one would be a > 3.6a anomaly (and yet such disks 
are found around <~ 2% of stars). However, given cur- 
rent limits in definitive detection of transient events, the 
derived incidence is probably a lower limit. 

One particularly interesting outcome of the model is 
the distribution of radii of the planctesimal belts. The 
model assumed a power law distribution of radii in the 
range 3-120 AU, and found from the distribution of radii 
of the sample detected at 24 and 70 /im, that this fol- 
lows N(r) oc r -°- 8 ± - 3 . it was found to be important to 
take the detection bias into account when interpreting 
the observed distribution. The origin of this distribu- 
tion is not considered in this model, which assumes that 
all stars have relatively narrow planetesimal belts with 
dr/r = 0.5. The inner hole in these belts has been in- 
ferred in other studies to be caused by the presence of 
inner planets (e.g., Roques et al. 1994; Wyatt et al. 
1999; Wilner et al. 2002; Wyatt 2003; Quillen 2006). 
If these are Kuipcr belt analogues (Wyatt et al. 2003), 
then this radius distribution would be indicative of the 
orbital radius of the outermost planet in its planetary 
system. However, these could also be asteroid belt ana- 
logues, i.e., planetesimal belts in the midst of a planetary 
system. If so the interpretation of this radius distribu- 
tion is less clear, and could, like the solar system for 
example, be indicative of the orbital radius of the most 
massive giant planet in these systems. In any case, this 
distribution provides an important and unique constraint 
on the outcome of planet formation models. 
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Another particularly powerful application of the model 
is that, regardless of the physical reason, it does explain 
the far-IR properties of A star disks and so it can also be 
used to make predictions for what surveys at other wave- 
lengths (or with different detection thresholds) will see. 
Such predictions are made much easier when the survey 
detection threshold is uniform, or at least can be read- 
ily characterized, and this study emphasizes the impor- 
tance of a characterizable threshold for future surveys. 
One such survey is the SCUBA-2 debris disks survey 
(Matthews ct al. 2007), which will include an unbiased 
sub-mm survey of the 100 nearest A stars to a uniform 
3cr sensitivity of 2 mJy at 850 /im. The model predicts 
that 17 of the 100 will be detected above 2 mJy, including 

5 above 10 mJy, and that 5 of the detected disks should 
be resolvable in sub-mm imaging on size scales > 6". 
However, the detected fraction could be higher than this, 
since it depends on the radius distribution, which these 
observations would set strong constraints on. It is also 
unknown how many A stars possess disks too cold to de- 
tect in the far-IR (Wyatt, Dent & Greaves 2003; Najita 

6 Williams 2005). This suggests that the SCUBA-2 sur- 
vey will be particularly fruitful, and the model presented 
here provides a framework which can be used to interpret 
the results of this survey, and those of future surveys. 



We are grateful for support provided by the Royal So- 
ciety (MCW), PPARC (RS, JSC) and SUPA (JSG). This 
work was partially supported by JPL/Caltech contract 
1255094 to the University of Arizona. 



Haisch, K. E., Lada, E. A., & Lada, C. J., 2001, ApJ, 553, L153 
Holland, W.S., et al., 1998, Nature, 392, 788 

Holmes E.K., Butner, H.M., Fajardo-Acosta, S.B., & Rebull, L.M., 

2003, AJ, 125, 3334 
Kalas, P., & Jewitt, D., 1995, AJ, 110, 794 

Kalas, P., Graham, J.R., & Clampin, M., 2005, Nature, 435, 1067 

Kenyon, S.J., & Bromley, B.C., 2002, AJ, 123, 1757 

Kenyon, S.J., & Bromley, B.C., 2004, AJ, 127, 513 

Krivov, A.V., Sremcevic, M., & Spahn, F. 2005, Icarus, 174, 105 

Lecar, M., Franklin, F.A., Holman, M.J., & Murray, N.J., 2001, 

ARA&A, 39, 581 
Liou, J.-C, & Zook, H.A., 1999, AJ, 118, 580 
Matthews, B., et al., 2007, PASP, submitted 

Moerchen, M.M., Telesco, CM., Packham, C, Kehoe, T.J. J., 2007, 

ApJ, 655, L109 
Moro-Martm, A., & Malhotra, R., 2003, AJ, 125, 2255 
Najita, J., & Williams, J. P., 2005, ApJ, 635, 625 
Natta, A., Grinin, V., & Mannings, V., 2000, in Protostars and 

Planets IV, eds. V. Mannings, A. P. Boss & S.S. Russell (Tuscon: 

Univ. Arizona Press) 
Nesvorny, D., Bottke, W.F., Lcvison, H.F., & Dones, L., 2003, ApJ, 

591, 486 

Quillen, A.C 2006, MNRAS, 372, L14 

Rhee, J.H., Song, I., Zuckerman, B., & McElwain, M., 2007, ApJ, 

in press, astro-ph/0609555 
Rieke, G.H., et al., 2005, ApJ, 620, 1010 

Roques, F., Scholl, H., Sicardy, B., & Smith, B.A. 1994, Icarus, 
108, 37 

Schneider C, et al., 2006, ApJ, 650, 414 

Sheret, I., Dent, W.R.F., & Wyatt, M.C, 2004, MNRAS, 348, 1282 
Song, I., Caillault, J.-P., Barrado y Navascues, D., & Stauffer, J.R., 

2001, ApJ, 546, 352 
Song, I., Zuckerman, B., Weinberger, A. J., & Bccklin, E.E., 2005, 

Nature, 436, 363 

Spanglcr, C, Sargent, A.I., Silverstone, M.D., Becklin, E.E., & 
Zuckerman, B., 2001, ApJ, 555, 932 



Steady-state evolution of debris disks around A stars 



19 



StaufTer, J.R., Hartmann, L.W., & Barrado y Navascues, D., 1995, 

ApJ, 454, 910 
Stern, S.A., 1996, AJ, 112, 1203 
Su, K.Y.L., et al., 2005, ApJ, 628, 487 
Su, K.Y.L., ct al., 2006, ApJ, 653, 675 
Telesco, CM., ct al., 2000, ApJ, 530, 329 
Telesco, CM., et al., 2005, Nature, 433, 133 

Weinberger A. J., Rich, R.M., Becklin, E.E., Zuckerman, B., & 

Matthews K., 2000, ApJ, 544, 937 
Williams, J. P., & Andrews, S.M., 2006, ApJ, 653, 1480 
Wilner, D.J., Holman, M.J., Kuchner, M.J., & Ho, P.T.P. 2002, 

ApJ, 569, L115 
Wyatt, M.C., 2003, ApJ, 598, 1321 
Wyatt, M.C., 2005, A&A, 433, 1007 



Wyatt, M.C, & Dent, W.R.F, 2002, MNRAS, 334, 589 
Wyatt, M.C, Dent, W.R.F., & Greaves, J.S., 2003, MNRAS, 342, 
876 

Wyatt, M.C, Holland, W.S., Greaves, J.S., & Dent, W.R.F., 2003, 

Earth, Moon, and Planets, 92, 423 
Wyatt, M.C, Dermott, S.F., Telesco, CM., Fisher, R.S., Grogan, 

K., Holmes, E.K., & Piha, R.K., 1999, ApJ, 527, 918 
Wyatt, M.C, Greaves, J.S., Dent, W.R.F. , & Coulson, I.M., 2005, 

ApJ, 620, 492 

Wyatt, M.C, Smith, R., Greaves, J.S., Beichman, C.A., Brydcn, 
C, Lisse, CM., 2007, ApJ, 658, 569 



